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We propose a Lie-algebraic duality approach to analyze non-equilibrium evolution of closed dy- 
namical systems and thermodynamics of interacting quantum lattice models (formulated in terms of 
Hubbard-Stratonovich dynamical systems). The first part of the paper utilizes a geometric Hilbert- 
space-invariant formulation of unitary time-evolution, where a quantum Hamiltonian is viewed as 
a trajectory in an abstract Lie algebra, while the sought-after evolution operator is a trajectory in 
a dynamic group, generated by the algebra via exponentiation. The evolution operator is uniquely 
determined by the time-dependent dual generators that satisfy a system of differential equations, 
dubbed here dual Schrodinger-Bloch equations, which represent a viable alternative to the conven- 
tional Schrodinger formulation. These dual Schrodinger-Bloch equations are derived and analyzed on 
a number of specific examples, ft is shown that deterministic dynamics of a closed classical dynam- 
ical system occurs as action of a symmetry group on a classical manifold and is driven by the same 
dual generators as in the corresponding quantum problem. This represents quantum-to-classical cor- 
respondence. In the second part of the paper, we further extend the Lie algebraic approach to a wide 
class of interacting many-particle lattice models. A generalized Hubbard-Stratonovich transform is 
proposed and it is used to show that the thermodynamic partition function of a generic many-body 
quantum lattice model can be expressed in terms of traces of single-particle evolution operators 
governed by the dynamic Hubbard-Stratonovich fields. The corresponding Hubbard-Stratonovich 
dynamical systems are generally non-unitary, which yields a number of notable complications, includ- 
ing break-down of the global exponential representation. Finally, we derive Hubbard-Stratonovich 
dynamical systems for the Bose-Hubbard model and a quantum spin model and use the Lie-algebraic 
approach to obtain new non-perturbative dual descriptions of these theories. 



I. INTRODUCTION 

The fundamental problem of quantum-to-classical cor- 
respondence is an old question that has existed since 
the birth of the quantum theory, and which appears 
frequently in essentially all fields of physics under var- 
ious covers, in particular, in the context of the path in- 
tegral approach to quantum mechanics of Feynman pQ, 
Glauber's coherent states [2] and their various generaliza- 
tions [3 5j, integrable systems [5J, classical and quantum 
chaos [7], AdS/CFT correspondence jS], and more gen- 
eral Langlands duality [5], etc. One particularly simple 
variation of this general problem is to connect classical 
motion of a particle or an ensemble of classical particles 
under the influence of external time-dependent fields to 
unitary evolution of the corresponding externally-driven 
quantum system. In many cases, including the text-book 
harmonic oscillator problem, the one-to-one correspon- 
dence between classical and quantum dynamics has been 
long known and well-established. However, the situa- 
tion remains less clear in some other cases, notably in 
dynamical systems where classical motion occurs on a 
manifold with constraints, which on the quantum side 
usually manifest themselves through topological terms in 
a quantum action or non-trivial commutation relations 
in the Hamiltonian formalism [TU]. The simplest such 
examples are an externally driven spin [11H17) or a two- 
photon system with time- varying parameters |18j . Vari- 
ous quasiclassical approaches [T9H22] have been proposed 
to analyze such systems, and the results of these analy- 



ses are often suggestive of classical and quantum unitary 
dynamics representing two sides of the same coin, but no 
simple general summary of such correspondence appears 
to exist to encompass all relevant cases on a systematic 
basis. 

Apart from the pure fundamental interest to the prob- 
lems of unitary quantum evolution and quantum-to- 
classical correspondence, the recent experimental ad- 
vances in the fields of quantum optics, atomic and molec- 
ular physics [25HTT] , and quantum control fET-pig] make 
these fundamental questions also practically important 
and experimentally relevant. Models that had existed 
before only in theorist's imagination can now be realized 
and studied in the laboratory, where unitary relaxation- 
less dynamics can be directly detected and analyzed. On 
the theory side however, the analytical tools to study 
such non-equilibrium dynamics are still being actively 
developed 0^53]. 

Motivated by these recent advances and our recent the- 
oretical results , we use the latter as well as results 
from the old [58H6T] and more recent [62H66] mathemat- 
ical physics literature and quantum control theory [51] 
to develop a general Lie-algebraic framework that de- 
scribes both unitary quantum dynamics and classical 
non-equilibrium evolution on an equal footing. We uti- 
lize the approach pioneered by Magnus [551 ES] and de- 
rive a set of differential equations, referred to below as 
dual Schrodinger-Bloch equations (DSBE), that make 
quantum-to-classical correspondence manifest. The ma- 
jor point emphasized throughout the paper is that the 
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unitary quantum evolution is largely decoupled from the 
Hilbert space and can be analyzed completely indepen- 
dently of it based on the minimal Lie algebraic structure 
provided by the underlying Hamiltonian. More specifi- 
cally, we first consider a Hamiltonian expressed in terms 
of operators that close under commutation into a finite- 
dimensional Lie algebra, A, which is specified by its 

d 

structure constants, \j a , J&] = i Yl f a b C ^c- The uni- 

c=l 

tary quantum evolution operator, U(t) lies in the Lie 
group, Ga — exp(iA), which is generated by the alge- 
bra via exponentiation and it is uniquely determined by 
a trajectory, & a (t), in the dynamical group as follows 



U(t) 



exp 



The classical equations of 



motion for the quantum-mechanical averages, M a (t) — 
(^(t)\J a \^{t)) follow from the structure constants of the 

d 

algebra, M b (t) = £ f ba c b a (t)M c (t). Classical motion 

a,c— 1 

generally takes place on a classical manifold with con- 
straints, A4, via action of a symmetry group, G(M), 
in the d-dimensional adjoint representation, M(t) = 

e~5-"* M(0), with f a representing the correspond- 
ing generators. It is demonstrated that this group action 
is determined by the very same dual generators, 
that appear in the quantum problem. This manifests 
quantum-to-classical correspondence. The dual gener- 
ators, $ a (i), satisfy a system of non-linear differential 
equations, the DSBE. New results derived here include a 
specific form of these dual Schrodinger-Bloch equations 
for a number of dynamical systems. 

The second part of the paper focuses on a class of 
more complicated Hamiltonians, bilinear in algebra gen- 
erators, that encompasses a wide variety of interacting 
many-body quantum systems of interest to both con- 
densed matter physics and cold atoms, which include 
various quantum lattice spin models and Hubbard-type 
models. We utilize the path integral formalism to for- 
mulate a generalized Hubbard-Stratonovich decomposi- 
tion [54] and to prove that thermodynamic properties 
of such many-body systems can be represented as a lin- 
ear combination of traces of quantum evolution operators 
corresponding to all possible realizations of the Hubbard- 
Stratonovich fields. This approach gives rise to Hubbard- 
Stratonovich quantum dynamical systems, which gener- 
ally involve non-unitary evolution and require a complex 
extension of the underlying Lie algebra. It is argued and 
demonstrated on an explicit example of the familiar Bosc- 
Hubbard model that in treating interacting many-body 
Hamiltonian, it is useful to Hubbard-Stratonovich decou- 
ple the terms associated with single-particle Lie algebra 
generators instead of focusing only on the interacting 
terms quartic in creation/annihilation operators in the 
Fock basis; i.e., the hopping terms, tf(i)b(j), and inter- 
action terms $ {i)b(i)b^ (j)b(j) represent the same level of 
complexity from this point of view and are to be decou- 
pled into linear combinations of b(i) 7 and b^(i)b(i), 



which together with the identity span the solvable four- 
dimensional harmonic oscillator algebra, f)4. We solve 
the corresponding Hubbard-Stratonovich dynamical sys- 
tem exactly and obtain a new non-perturbative repre- 
sentation for the partition function of the Bosc-Hubbard 
model. We also analyze a frustrated quantum spin model 
and, using the dual Schrodinger-Bloch equation, derive a 
dual representation of the theory in terms of group gen- 
erators. It is argued that such dual descriptions may 
be helpful in circumventing the sign problem in certain 
cases. 

Our paper is structured as follows: In Sec. [TTJ we, 
starting from the standard Schrodinger equation in a 
Hilbert space, derive a Hilbert-space-invariant form of 
the Schrodinger equation for quantum evolution that 
relies only on the primitive Lie-algebraic structure of 
the dynamical system. In Sec. |III[ we utilize the 
global exponentiation conjecture (that quantum evolu- 
tion can be globally exponentiated from the algebra, 
U(t) = exp [— i3>(i) • J] ) and derive a general form of the 
dual Schrodinger-Bloch equation for the Hilbert-space- 
invariant generators, Sec. |IV| provides a number 
of examples of the usage of the dual Schrodinger-Bloch 



equations: Sec. IV A derives and solves DSBE for the 
harmonic oscillator algebra, Sec. IV B| derives DSBE for 
the su(2) spin algebra, Sec. IV C establishes a relation 
between the latter two dynamical systems, Sec. [IVD] dis- 
cusses the two-photon algebra, and Sec. IV E formulates 
a generic five-step procedure of reverse-engineering ex- 
act solutions using DSBE in the higher-rank su(N) dy- 
namical systems. Sec . [V] establishes quantum-to-classical 



correspondence and Sec. |VB| illustrates the meaning of 
this correspondence on the simple example of linearly- 
driven Jaynes-Cummings model that is solved exactly. 



Sec. VI extends the Lie-alegraic ideas to interacting 
many-particle systems: a generic Lie-algebraic lattice 
model is introduced in Sec. IVI Al and Sec. IVI Bl derives a 
generalized Hubbard-Stratonovich transform and formu- 
lates a generic Hubbard-Stratonovich dynamical system. 



Sec. VII includes two examples of this Lie-algebraic ap- 
proach to interacting quantum models: In Sec. |VII A| we 
solve exactly the Hubbard-Stratonovich dynamical sys- 
tem for the harmonic oscillator algebra to derive an ex- 
act, previously unknown, representation of the partition 
function of the Bose-Hubbard model. Sec. IVIIBI formu- 
lates, on the example of a frustrated quantum spin model, 
a dual path-integral approach to interacting spin systems 
where the integration over Hubbard-Stratonovich fields is 
replaced with path-integration over the dual generators 
in the dynamical group. Sec. |VIII| provides a summary 
and discussion of many open questions. Appendices [A] 
and [5] are logically connected to the first "single-particle" 
part of the paper, namely they relate to Sees. |III| and |IV| 
and Sec. |VJ correspondingly: Appendix [X] summarizes 
general algebraic properties and classification of possible 
unitary quantum dynamics and appendix [5] reviews the 
construction of the classical manifold that plays a central 
role in the discussion of quantum-to-classical correspon- 
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dence. 



II. FROM SCHRODINGER PICTURE TO A 
HILBERT-SPACE-INVARIANT FORMULATION 

A. Dynamical Lie algebra 

Consider the non-equilibrium Schrodinger equation 



id t Mt))=H L (t)\*(t)) 
l*(0)) = l^o), 



(1) 



where "Hl(^) is an L x L Hermitian Hamiltonian matrix 
acting on the wave-functions, in an L-dimensional 

complex Hilbert space 'Hil(L) spanned by the linear com- 

L 

binations, — c k\k) , where \k) G Hil(L) are or- 
fc=i 

thonormal basis vectors in the Hilbert space and cjEC. 
\ipo) = X] c| \k) G Hil(L) is a normalized initial state. 

k=l 

Note that the normalized wave- functions lie on the hyper- 
surface of the (2L — l)-dimensional sphere, S 2i_1 . In this 
Section, we assume that L < oo but the key results will 
be generalized to a specific class of infinite-dimensional 



representations in appendix A 4 



Any such dynamic Hamiltonian can be written in 
the form, 7^i(t) = P % (t)®i, where Oi are Hcrmi- 

i=l 

tian linearly-independent matrices and P l {t) are time- 
dependent parameters. Also, assume, without loss of 
generality, that TrC); = 0. If this is not so from the 
outset, we can always redefine our operators as fol- 
lows 6i -» (6i- ^TrdX where I L is the L x L 

identity matrix, and trivially solve the Abelian part 
of the Schrodinger equation that will give rise to a 
pure phase dynamics of the wave-function as follows, 



!*(*)> 



exp 



| *(*)). All oper- 



i ^TrO i J(3 i (s)ds 
»=i o J 
ators in the remaining part of the Hamiltonian will be 

traceless and therefore our assumption does not reduce 
the practical generality of the analysis. This implies that 
these traceless operators, Oi, certainly belong to an L- 
dimensional defining representation of the special linear 
algebra st(L,C) and the unitarity constraint further re- 
quires the operators to belong to the L-dimensional uni- 
tary representation of su(L). 

However, the quantum dynamical system may gen- 
erate a much smaller Lie algebra than the (I? — 1)- 
dimensional su(L) (recall that L is the dimensionality 
of the Hilbert space). In fact, we argue below that this 
algebra, A, often has little-to-no relation to a particular 
Hilbert space. Consider now all possible matrix com- 



mutators Oij = —i 



Ou O, 



and multiple commutators, 



, etc. that are guaranteed to eventu- 



(or more precisely into its L-dimensional representation, 
T L [A]). Note that m < d = dirndl < (L 2 - 1). 

Choose a basis in [A], i.e., a set of d linearly- 
independent (L x L) Hermitian traceless matrices, J a £ 
Tl [A], which themselves are linear combinations of ma- 
trices Oi in the original Hamiltonian and their multiple 
commutators. These generators will satisfy a specific set 
of commutation relations, 



(2) 



Jai Jb\ — / fab C J< 



where the structure constants are guaranteed to be real, 
f ab c G K. Finally, we can re- write the matrix Hamilto- 
nian in terms of the generators as follows 



(3) 



Eqs. Q and ^ summarize the fact that an arbitrary 
time-dependent Schrodinger equation formulated in a 
finite- dimensional Hilbert space can be reduced to the 



Lie-algebraic form (see also, appendix A 4 for a discussion 



of certain theories in infinite-dimensional Hilbert spaces) 



B. Exponential image and dynamical group 

Lie algebras usually appear as a means to study Lie 
groups, and are often considered by-products of the lat- 
ter. In the context of non-equilibrium quantum mechan- 
ics however, a reverse view is much more useful |51j . Be- 
low, we treat the Lie algebra, A, arising from the under- 
lying Hamiltonian as a primary structure, while the Lie 
group it generates, Ga — exp(iA), as secondary. A par- 
ticular Hilbert space will be a tertiary structure of minor 
or no importance. To proceed with this program, we use 
structure constants in Eq. ^ to define an abstract Lie 
algebra, A, spanned by the generators, J a , 



[Ja, Jb] = ' Ad J a Jb = i.fab C Jc 



(4) 



ally close into a finite-dimensional Lie algebra, A C su(L) 



where here and below a sum over repeating indices is as- 
sumed, unless stated otherwise. Here we put inverse hats 
on top of the abstract generators, J a , of the abstract al- 
gebra, A = span {Ji, . . . , Jd} , to distinguish them from 
their specific matrix representation arising from the con- 
ventional Schrodinger formulation in the Hilbert space, 
'Hil(L). We note that the word, "abstract," here implies 
that J's are not matrix objects, nor are they operators 
acting on a wave-function, but rather vectors in a lin- 
ear space equipped with the commutator that makes it 
into a Lie algebra. It also implies that the products, 
JaJb, etc- do not make sense in the abstract Lie 
algebra, but only commutators do. The products are 
of course well-defined given a representation but nor- 
mally have little relation to the underlying dynamical 
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system. E.g., the square of an su(2)-generator from its 
two-dimensional representation (i.e., half- integer spin de- 
scribed by the Pauli matrices) is a^. y z /4 = 7 2 /4 and is 
proportional to an identity matrix that commutes with 
any other observable. On the other hand, the square of a 
generator from any higher-dimensional representation of 
su(2) (e.g., spin-1 particle described by 3 x 3 matrices) is 
non-trivial and as such may in principle have non-trivial 
dynamics. In what follows, we mostly concentrate on 
Lie-algebraic invariants, which are defined as objects or 
statements, whose definition or proof relies only on the 
commutators of the generators in the dynamic Lie alge- 
bra and their linear combinations and does not involve 
any other algebraic structures. 

However, a particular faithful representation can be 
used as a tool to prove certain Lie-algebraic-invariant 
statements and identities and demonstrate that they are 
in fact independent of the representation itself. The 
Hadamard lemma (think, unitary rotation of an operator 
into the Heisenberg or interaction picture), 

and the Baker-Campbell-Hausdorff (BCH) relation for 
e tX e tY = e lZ (think, terms in a time-ordered exponen- 
tial), 

l 

Z = X + [ dsij) (e iad * e ls ad ^ ) Y, (6) 



stand out as perhaps the most important among such in- 
variant identities [in Eqs. (5]) and ([6]), we assumed that 
X ,Y £ A and used the logarithmic derivative of the 
Gamma-function, defined via its series, e.g., ip(e z ) = 

oo 

^fz n , with B n being the Bernoulli numbers]. 

n=0 

While the usefulness of Eq. Q for practical calcula- 
tions appears to be limited, it does provide a nice il- 
lustration of the invariance of the BCH-product. Let 
us recall however that in the case of an arbitrary finite- 
dimensional Lie algebra, the BCH identity is convergent 
in a close vicinity of the origin, G A, but Eq. ^ does 
not always converge globally 67J, with the algebra as sim- 
ple as sl(2, C) being the classic counter-example. Another 
manifestation of this problem is that there is generally 
no guarantee that in a particular matrix representation 
a product of two matrix exponentials of algebra elements 
can be written as a single matrix exponential of an alge- 
bra element. This "complication" is also closely tied to 
the fact that a Lie algebra is guaranteed to exponentiate 
into its Lie group only locally, but does not always do so 
globally. 

One can argue however that in the context of quan- 
tum dynamical system ([I]) , the dynamic algebra does ex- 
ponentiate onto a Lie group. Here we recall that the 
covering problem of global exponentiation of an arbi- 
trary finite-dimensional Lie algebra arises due to non- 
compact generators that do not necessarily form a sub- 
group [68 . However, such generators are prohibited 



in the quantum dynamical system described by Hermi- 
tian finite-dimensional matrices. The very existence of a 
finite-dimensional Hilbert space guarantees that the set 
of matrix exponentials and their products, form a subset 
of SU{L) group, which is indeed compact. This sug- 
gests that the entire dynamic group is a compact man- 
ifold that can be globally generated by an exponential 
map. An argument here is that any two points of this 
compact Riemannian manifold are connected by a min- 
imizing geodesies and the exponentials from the algebra 
generate such geodesies. Therefore, it seems that we can 
safely define a compact abstract group, generated by the 
abstract dynamic algebra via exponentiation that we de- 
note as 



G 



A 



JA 



(7) 



We emphasize that while this dynamical symmetry group 
may be a subgroup of SU(L), it generally has little to do 
with a particular Hilbert space, which we have used here 
only to argue in favor of the existence of a global ex- 
ponential map onto the Lie group, Gj±. We will refer 
to the assumption about the existence of this map as to 
"global exponentiation conjecture" (we call it a "conjec- 
ture," because the mathematical rigor of the above ar- 
guments is not entirely satisfactory and the background 
of the author is insufficient to be able to formulate a 
general statement in the form of a rigorous theorem). 
Finally, let us stress that this conjecture does not as- 
sume that the specific form of BCH formula ^ and/or 
Magnus series [SSI US] are globally convergent, but it 
does imply, in particular, the following: Given a faithful 
matrix representation, the matrix product of two expo- 
nentials of algebra elements gives rise to an exponential 
of another algebra element, i.e., if A", Y £ Tl[A], then, 
-i\n(e l *e lt ) £ T L [A]. 



C. Geometric Alternative to Unitary Evolution in 
a Hilbert space 

Returning to our quantum dynamical system ([I]), we 
observe however, that since the familiar Schrodinger 
equation is written in terms of a wave-function, it is man- 
ifestly representation-specific. Let us now cast it into a 
representation-invariant form. In the framework of the 
Schrodinger formulation (JlJ, our goal is to find a norm- 
preserving solution governed by a unitary Lx L evolution 
matrix as follows 

L 

mt)) = u L (t)\ih)= E c* p (*)4 0) i*>- ( g ) 

k,p=l 

We can now substitute this solution in the Schrodinger 
equation to arrive to 



id t U L {t)=H L (t)U L {t) 

Ul(0) = h. 



(9) 
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fonian 



■lives 



Compact Lie group, G=exp(iA) 
. (quantum evolution occurs 



'(O 



is a trajectory in the abstract Lie algebra, A, while the 
evolution operator is a trajectory in the abstract Lie 
group, U(t) C G_4 = exp[£A]. To solve the quantum- 
dynamical system ( 10 1 implies to find the latter from the 
former. Fig. [I] illustrates this geometric view of quantum 
dynamics. 



III. DUAL SCHRODINGER-BLOCH EQUATION 




FIG. 1: Cartoon illustration of the geometrical view of quan- 
tum dynamics: The Hamiltonian is a trajectory in a Lie al- 
gebra, A, which is a finite-dimensional vector space. This 
trajectory gives rise to unitary evolution, which is a trajec- 
tory in the Lie group, Ga = <^ A , generated by the algebra via 
exponentiation. Topologically, this group is a compact con- 
nected manifold. As emphasized in the text, neither A nor 
Ga are associated with any particular Hilbert space. 



As evident from (|9| , the quantum evolution operator is 
decoupled from the initial condition, which is well-known. 
What is also known [55] , but perhaps less widely appreci- 
ated, is that the unitary evolution is also decoupled from 
the Hilbert space and is representation-invariant. To see 
this, multiply H by U£ _1 (t) = Jjl(t) from the right to 

find id t V 'i,(t)Ujjt) = ti^t)- Based on the global expo- 
nentiation conjecture of the previous section, we conclude 
that since at t = the identity, Ul(0) = II, certainly 
belongs an L-dimensional representation of the group, 
Tl [Ga], it w iU stay there at all times. The same conjec- 
ture also implies that the solution can be sought in the 



form 58 , Ul (t) = exp 



where <!> a (t) are real 



-i$ a {t)J a 

functions of time [we will also use the vector notation, 
* = ...,&), J = (j x , J 2 ,..., J d y and a dot- 

product <&(t) ■ J = $ a (t)J a ]. Now we observe that the 
product, idtUL(t)U^(t), is representation-invariant, be- 
cause it can be calculated using the BCH-identity only. 
Therefore, one can drop the representation indices and 
write the Schrddinger equation in an "abstract form:" 

id t U(t)U- l (t) = i \d t e-^^- J } e **(*M = b a (t)J a 



U(t) G G A and U(0) = 1 



jo 



(10) 

Note that eventhough, we have "derived" this equation 
from the textbook Schrodinger equation ([I]), Eq. ( 10 ) has 
no "recollection" about the Hilbert space and can be an- 
alyzed completely independently of it. In fact, we might 
have as well started from Eq. ( 10 1, which does not involve 



the notion of a wave-function at all. In this approach, the 
Hamiltonian 



The evolution operator in Eq. ( 10 1 is a trajectory in 



the Lie-group, generated by the Lie-algebra of the driving 
Hamiltonian, and therefore has the form 

U(t) = exp [-i$ a {t)J a ] and J7 _1 (i) = exp [+i$ Q (i)J a ] . 

(12) 

We note that since the vector, <&(£), defines the group, 
there appears a problem of its proper parameterization 



that we will discuss at a later stage (see, Sec. IVF2| 



However, we shall not dwell on this subtlety for the mo- 
ment and treat the corresponding differential equation 
( 10 1 rather freely, keeping in mind that there exists a 



finite-dimensional matrix representation that always can 
be used to justify the corresponding algebraic manipula- 
tions. 



Eqs. ( 10 ) and ( 12 ) lead to the following structure, 

i lim < - 

-s-o I e 



H(t) = ilimlMe-^'W+^Wl^e^'W^-e* )}, 

which had been considered before by many, includ- 
ing notably Feynman (see, e.g., Ref. [59] for a re- 
lated discussion). Using the Hadamard lemma ([5]) and 

the related operator/matrix identity, §^e x+eY 
Jq 1 dse^-^XYe 3 * , we obtain 



<S> a {t) J dse- ls * (t)adj fc -6 Q (i) 
" 



Ja =0, 



(13) 



which is to be supplemented by the initial condition, 
= 0. We see that Eq. (13 1 involves nothing but 



commutators and therefore is representation invariant in 
the sense that the vector function, i&(t), determines evo- 
lution operators in all Hilbert spaces. 

The operators, (ad # .j), naturally give rise to a d- 
dimensional matrix representation of A, which is known 
as the adjoint representation of the algebra, and it is 
uniquely determined by the structure constants. The 
representation is constructed in a standard way by as- 
signing to each vector the following (d x d)-matrix, 



i(ad # .j) J b = * ■ f b c J c 



(14) 



H(t) = b(t) ■ J 



(11) 



The matrices, fi,---fd correspond to the adjoint rep- 
resentation of algebra generators Ji, . . . ,Jd- Notation- 
wise, here and below, the hat on top of the /-matrices 
signifies the fact that they are indeed matrices (rather 
than abstract objects) and the bold fonts correspond to 



G 



Note 



d-dimensional vectors, e.g., / = (fi, fo, ■ ■ ■ , f 

also that if we go back to the matrix elements, f ab c we re- 
cover the structure constants of the algebra by construc- 
tion. Therefore, the adjoint representation is basically 
a set of d X d matrices built directly from the structure 
constants. Using these matrix notations, we can write 
Eq. ( 13 ) in the following compact form 



in universal simple algebras (for practical purposes, un- 
derstanding driven su(L) systems should be sufficient) 
and dynamics in non-universal solvable algebras (that are 
argued to be comparatively trivial and always exactly 
solvable via factorization). In this section, we demon- 
strate the application of the DSBE for both solvable (har- 
monic oscillator) and simple [su(L)] algebras. 



dse'* w -'*(t) = b(i), *(0) = O 



(15) 



This is the main equation used throughout the paper. We 
will refer to it as the dual Schrddinger-Bloch equation 
(DSBE), and use it for actual calculations of quantum 
and classical dynamics. Let us emphasize two important 
properties of the DSBE, which will be elaborated upon 
in the following sections: (i) Eq. ( 15 ) has no "recollec- 
tion" of the original Hilbert space and in its explicit form 
does not contain the Planck constant; (ii) DSBE is purely 
real, i.e., / are real- valued matrices and the sought-after 
functions are real as well. 

A word of caution is in order: Since any matrix, in- 
cluding (3> • /) in Eq. (15), commutes with itself, it is 



tempting to calculate the integral over s in (151 and 



rewrite it terms of the inverse matrix and its exponen- 
tial, c.f., Refs. [B21|53]. However, one should avoid doing 
so and exercise care at this point, because the adjoint- 
representation matrix does not necessarily have an in- 
verse, i.e., = oo. We emphasize that this is not a 
minor point, but an essential complication, which should 
be kept in mind when analyzing the DSBE. A correct 
approach is to calculate the matrix exponential, e s *^^'^, 



first, and evaluate the integral over s in Eq. (151 only 



afterwards (reversing the order of these operations may 
lead to errors in actual calculations). 

Concluding this section, let us also recall that the 
adjoint representation defines the Cartan-Killing inner 
product in the algebra, 



(*, 0)=Tr {(*■/) (©'/)} 



(16) 



which is used as a tool to classify the algebras 68 . This 
Cartan-Killing product is non-singular for semisimple al- 
gebras and in this case the Cartan-Killing product of ba- 
sis vectors, can be identified with the metric that gives 
rise to geometric structures. 



A. DSBE for the Weyl-Heisenberg algebra, f)4 

Consider the four-dimensional Weyl-Heisenberg alge- 
bra f)4 = span {l, d, d^, d^d} = span {l, x,p, d^d}, asso- 
ciated with a driven harmonic oscillator with the Hamil- 
tonian 

U^Jt) =Lj(t)tfa + a*(t)a* + a(t)a + b°(t)l, (17) 

where the operators a = aT and a) = d + satisfy 
the canonic commutation relations, [d, aH = 1 and 
\a*a, d^ = id 1 * 1 and x = (a + <v)/\2 and p = 
i(d^ — a)j\pl satisfy [x, p] = il. Want to find the so- 



lution to Eq. (I10j> corresponding to (17 1. To illustrate an 



application of the general factorization scheme outlined 
in appendix |A 2[ let us decompose the solvable algebra 
1)4 into an Abelian and nilpotent components as follows 

f) 4 = span {d^d} + span {I, x,p} = u(l) + f) 3 , (18) 

where the latter, is a nilpotent subalgebra. Solve the 
other u(l)-component by introducing the evolution ma- 



trix U u (i)(t) = exp 



-iata J uj{t')dt' 



. The full evolution 



operator is given by the product U\ ]i = t7 u (i)(i)t/(, 3 (i), 
where the last factor is sought in the form 

Ut, 3 (t) =exp{-i [$°(t)l + $ 1 (t)i + $ 2 (i)p]} (19) 

and is to satisfy the Schrodinger equation with the ro- 
tated Hamiltonian. 



H t)3 (t)=b°(t)i + b 1 (t)x + b 2 {t)p, 



(20) 



i f uj{t )dt 

where 6 x (t) + ib 2 (t) = ^a(t)e (recall that 

exp [— iJ7ad a ta] a = a ). Using the regular rep- 
resentation of f)3, we can present the operator (14) in the 



DSBE ( 15 ) as follows: 



IV. EXAMPLES OF THE DUAL 
SCHRODINGER-BLOCH EQUATION 

As argued in appendix [X| the variety of all possi- 
ble quantum dynamical systems formulated in a finite- 
dimensional Hilbert space and certain theories formu- 
lated in an infinite-dimensional Hilbert space (see ap- 



pendix A 4), reduce to an analysis of dynamical systems 





F 64 [*(t)] = #(t)-/ = ( $ 2 (U 



(21) 



Due to the nilpotent structure of the algebra, the matrix 
(21) exponentiates easily, leaving behind only two non- 



vanishing terms in the series, and yields f dse^^'f 

o 
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( 1 °\ 

$ 2 (£)/2 10. The dual Schrodinger-Bloch equa- 

\ -$ x (t)/2 1 / 
tion therefore reduces to 



$ x ' 2 (i) = b^ z {t) 



(22) 



and 



$°(t) + - [^(^(t) - $ 2 (t)$ x (t)] = 6°(t). (23) 
Hence, we find the exact quantum evolution as follows: 



, , ifu(t")dt" 



+ i$2 W = 77f / ^'«(*')e l ° 



(24) 



and 



t t' 



1 



= / *'6°(t') - ^ / di' / dt" \a(t')\ \a(t")\ 



x sin i^Jdr [w(t) + - 0(t")] J , (25) 

where a(t) = |a(f)| e^ ( *>. 

Finally, note that the phase space associated with the 
harmonic oscillator algebra is not compact, which seems 
to be in conflict with the arguments of Sees. |II A| and 
II C (see also Fig. [I]) , which argued that the dynamical 
group must be compact. The resolution of the paradox 
is in that the harmonic oscillator algebra does not have 
a finite-dimensional unitary representation and hence no 
finite-dimensional Hilbert space exists, which was a key 



assumption of Sec. II A On the other hand, the expo- 
nential representation still works perfectly well (see also, 
Sec. IV C and appendix A 4 1. 



B. DSBE for the su(2)-algebra 

Now consider the simple algebra, su(2), defined via the 
commutation relations 



(26) 



where the Greek indices run over 1,2,3 = x,y,z and 
here and below in this section we make no distinction 
between the covariant and contravariant indices. We are 
interested in spin dynamics governed by the Hamiltonian, 



#«u(2)(*) =b(t)- J. 



(27) 



We follow the general prescription outlined above and 
also in Ref. [57], where the su(2) case was considered 
in the context of two-level system dynamics, and seek 



the solution in the form U, 



0U(2) 



exp • J I . The 



adjoint operator in the DSBE (15), (— iad^^^.j 
regular matrix representation as follows 



F 



SU (2) 



[*(*)] 





f 


-*» 




* • / = 1 



















has a 



(28) 



and describes a generator of S'0(3)-rotations. The ex- 
ponential of this operator is therefore a finite rotation 
around the axis n = 3>/ < I > , and the exponential from 



Eq. ( 15 ) reads 



e s*-/ _ cos s< j) _|_ 7j Q 7 1(9 [1 _ C os s$] — e Q ^ 7 n 7 sin s$. 

(29) 

We now present the time-derivative as follows $ = 
<I>n + $n, and keeping in mind that h ■ n = 0, obtain the 
DSBE for the su(2) algebra 



$ n + sin <i> n + (1 - cos$) [nxn]= b(i). 



(30) 



This equation is not formally solvable for a generic time- 
dependence of the driving field. However, it could be 
used to prove certain exact relations [such as, e.g., $(i) = 
J Q drb(r) • n(r)] and generate exactly solvable models 
for two-level-system dynamics |llH17j by inverting the 
trajectories in the SU(2) group back onto the algebra. 



C. Contracting the DSBE from the u(2) algebra 
into Weyl-Heisenberg 

The harmonic oscillator algebra and spin algebra as 
well as their dual Schrodinger-Bloch equations discussed 
above appear to be completely different: driven dynam- 
ics in f)4 is solvable, while that in su(2) is not exactly- 
solvable; the former algebra is associated with the non- 
compact Euclidean classical phase space (x,p) € M 2 , 
while the latter with the magnetization dynamics on the 
Bloch sphere, M E S 2 . However, despite these drastic 
distinctions, the two algebras are known to be closely 
connected and this connection beautifully manifests it- 
self via Lie-algebraic contraction. We recap the main 
idea of this particular contraction [55] , because such pro- 
cedures generally provide an effective means to reduce 
Schrodinger-Bloch duals of an algebra to the correspond- 
ing equations in more primitive algebras that descend 
from it via contraction. 

Consider the four-dimensional u(2) algebra that is ob- 
tained from su(2) by adding to it a one-dimensional 
Abelian component, 



u(2) = span {Jo, J} = u(l) +su(2), 



(31) 



where [Jq, J a ] = 0. Define a parameter-dependent lin- 
ear transform of the u(2) generators: Jo — Jq, Ji,i{t) — 
eJ y ,x, and 1/3 (e) = J z + e Jo- The non-trivial com- 
mutation relations in terms of the new operators are 
[Js, 1/1,2] = =RJ2,i and [Ji, J 2 ] = -ze 2 J3 + ij . In 



the limit of e — > 0, these commutation relations become 
identical to those in the Weyl-Heisenberg algebra, f)4. 
I.e., 

limu e (2) ee lim span {j Q ,Jx{e),J 2 {e),J 3 {z)} 

— span {l,x,p, d^d} = f)4. (32) 



Note that the DSBE for u(2) follow from Eq. (30) by 



adding to it a trivial equation for the Abelian u(l) compo- 
nent. Let us however not solve it but instead re- write the 
"trajectory" in the algebra in terms of the new operators 
<J> a (i)J a = 9 a (i)Ja, which leads to 6° = <J>° - e" 2 * 2 , 

6 1 = $V e j O 1 = ® v / e > and e3 = $3 - The corre- 
sponding substitutions lead to the evolution operator, 



and new set of equations 



>'u e (2) = exp -i £ G a (t)J a 

a=0 

for 6°(i) that in the e — > limit crossover to the DSBE 
for fj4. 



D. Six-dimensional two-photon algebra, fie 



of the standard model. To calculate directly the expo- 



nentials of 8-dimensional matrices in the DSBE ( 15 ) is a 



cumbersome exercise and we proceed differently in a way 
generalizable to other higher-rank algebras. 

First, recall that there exists a three-dimensional 
(defining) representation of su(3) that can be constructed 



as follows. Define, t v 



~(l,r, 



with 1 < I < m < 3 



and a = x,y,z = 1,2,3 to be three-dimensional ma- 
trices such that the only possibly non-zero elements are 
(l, I), (/, to), (to, I), and (m, to) [here the notation implies, 
(row, column)] and the corresponding 2x2 matrices are 
the familiar Pauli matrices, a a . The standard Gell-Mann 
matrices, introduced in relation to the famous "eight- 
fold way," are expressed via r's as follows: Ai = f(x,2)jx) 



A 2 = 



r (l,2)-y, ^3 



A, = 



r (l,2);z, M 



A 4 = 



r (l,3);x> A 5 



A, = 



r (l,3);y; 



A6 = T(2,3);x, ^7 = T(2,3);j/j A 8 = [t(2,3);z + ^(l,3);z] ■ 

Note that only two among these Gell-Mann matrices are 
diagonal, A3 and As, and they correspond to the three- 
dimensional representation of mutually commuting Car- 
tan operators. We are interested in studying su(3) dy- 
namics driven by the time-dependent Hamiltonian, 



The two-photon operator algebra, f)6 = 
span {l, a, a\ a) a, d 2 , (d^) 2 }, is discussed in detail 
in the review by Gimlore et al. [18] Hence, we just reiter- 
ate the main results in a slightly different form relevant 
to the preceding discussion. \)§ is not semisimple, but in 
accordance with the general rule, it can be decomposed 
into a simple component and a nilpotent component as 
follows 



span \ - [{a)) 2 



span {l, x,p} = su(2) + f) 3 



(33) 



itl 



1. 



Using the canonic commutation relations, [d, 
we can verify that the first three generators in Eq. (|33|) 



give way to the commutation relation s id entical to those 
of the "usual" spin J x , y ,z operators (|26|. We can also 



explicitly verify that [su(2bl f) 3 ] C f)3 in (33 1. Hence, we 
can utilize the DSBE Eq. ( 30 ) to determine the dynamics 



of the first component and the factorization trick ( A3 ) 



or/and the DSBE for the nilpotent component (22 23) 
to determine the total evolution matrix, 



U fl6 (t) = U M2) (t)U t)3 (t)- (34) 

Reverse-engineering exact dynamics in 
higher-rank groups; SU(N) examples 

1. Dynamic "eightfold way" 



E. 



The next level of complexity after su(2) is provided 
by rank-2 algebras, such as su(3). We first focus on an 
su(3) dynamical system [51] associated with the latter 8- 
dimensional algebra, which is well-studied in the context 



T 3 [H su(3) (t)] =J2b a (t)\ a , 



(35) 



a=l 



where Ai 2,... 8 are the eight three-dimensional Gell-Mann 
matrices. Now consider the operator in the exponential 



(Va q ). 



of the DSBE (|15|) in this representation, i.e., the operator 

Va„ 



There exists a unitary time-dependent SU(3) 
rotation that diagonalizes this matrix into a combination 
of Cartan generators 

$ a (t)U 3 (t)\ a Ul(t) = wi(*)Aa + ^ 2 (<)A 8 . (36) 

Now, we can use the Bulgac-Kusnezov construction [70] 
or alternatively seek the unitary rotation as a product of 
three SU{2) rotations 

U 3 (t)= 11 expHx^-Ti,], (37) 

"=(*ij>0 

where we can take x"{t) to be "xy- vectors," i.e., x"(t) • 
t v = x v ' x {t)^v,x + X V]V {t)"^u;y Note that per elementary 
properties of the Pauli matrices, we can write 



Us(t)= II [^ocosl^WI-^n.^-f^sinlx^t)! 

v=(i,j>i) 

(38) 

where n u {t) = X u {t)/ \x v {t)\ i s a two-dimensional unit 
vector. We emphasize that unlike Eq. (37), Eq. (38) 



is not representation-invariant. However it is useful for 
actual calculations and can be utilized to relate $ a (t) 
and ujiflit) in Eq. (36), and these relations are indeed 



representation-invariant . 

Now, present the unitary rotation in the eight- 
dimensional adjoint representation using the invariant 



9 



form (37) and write the DSBE as follows 
1 



ul(t)U 8 (t) / dse s * a W?*tf s (t)U 8 (tWt) = b(t), (39) 



where f a are the known su(3)-generators in the eight- 
dimensional representation. This rotation brings the 
form <j> a (i)/ a to a diagonal form and we have 

l 

Ul(t) I dse s ^ lit)f3+ul2{t) ^U s (t)<b(t) = b(t). (40) 



Recall that f a=3 = diag {2, -2, 0, -1, 1, 1, -1, 0} and 
fa=s = \/3diag {0,0,0,1,-1,1,-1,0}. Therefore, 
the matrix exp onentiation of the diagonal matri- 
ces in Eq. (40) is simple and yields for T>(t) = 



/ dse-^M'^+^M/s] the fon ow i n g explicit expression 
o 

t>(t) = di&g { w[2uji], w* [2loi], l,w[VEui2 - wi], 
w*[V3uj 2 + Wl], 1}, 

where 

-2isz -iz Smz 



(41) 



w(2z) = I dse 



(42) 



Note that lim w(z) = 1. The third and eight matrix el- 

cmcnts that are identically e qua l to one in the diagonal 
matrix function T>(t) in Eq. (41) correspond to the zero 

eigenvalues of the generic matrix, ^<£> a / a ^, which there- 
fore does not have an inverse. It is a generic feature of 
the adjoint representation, which descends from the sim- 
ple fact that all elements in a Lie algebra commute with 
themselves. Finally, the Schrodinger-Bloch dual for su(3) 
can be written as 



Ul(t)V(t)U s (t)f>(t) = b(t), 



(43) 



where Us(t) is the eight-dimensional representation of 
the unitary SU(3) rotation (37) that diagonalizes the 

rotating three-dimensional "Hamiltonian, 
which has three real time-dependent eigenv alue s, 
[wa(t)/V3±wi(*)] and [-2u} 2 {t)/V^\ [see, Eqjje)], 



and the diagonal matrix D(t) is given by Eqs. (41) and 
( |42| ). Note that the main difficulty in deriving the dual 
dynamical system (43) is to find the frequencies, wi^i), 



and the rotation matrix, Us(t), with the latter to be ex- 
pressed in terms of the representation-invariant gener- 
ators. Technically, deriving these quantities reduces to 
solving a generic cubic algebraic equation, which is cum- 
bersome but should always be practically doable. This 
observation also suggests that an explicit derivation of 



the DSBE for an arbitrary simple algebra, A, is possible 
if the dimensionality of the minimal faithful representa- 
tion is i m i n < 4 and connects this problem nicely to the 
Galois theory. 

However, even if an explicit representation invariant 
dynamic su(3) system is written down, to solve it [i.e., 
to find eight dynamic generators, $ a (t)] may be a hope- 
less goal, as it generally amounts to the gargantuan task 
of solving a system of eight coupled differential equa- 
tions. On the other hand, just like in the su(2) case, 
the inverse problem of determining driving fields that 
give rise to certain non-linear evolution is much simpler. 
Eq. ( 43 ) should then be viewed as a result for the eight- 
dimensional "su(3) Zeeman magnetic field," b(t). 



2. Generalizing the results to su(N) 

The above discussion is not specific to su(3), but gen- 
eralizes to an arbitrary simple algebra and in particular 
to su(N). su(JV) is the most symmetric among all al- 
gebras and can be viewed as N(N — l)/2 su(2)'s tied 
together. The TV-dimensional defining representation of 
su(N) can be constructed in much the same way as in 
su(3), by introducing Pauli matrices, f„. Q with v = 
such that 1 < i < j < N and a — x,y, z. However, only 
(N — 1) linearly-independent diagonal generators (Car- 
tan operators, hi) exist, which is exactly the rank of the 
algebra, EHank su(N) = N — 1 = r. Therefore, the entire 
algebra is spanned by N(N — l)/2 Pauli matrices of the 
x and y-type (or equivalently, by raising/lowering oper- 
ators, T v - t ±) and the (N — 1) Cartan operators, which of 
course correctly reproduces the dimension of the algebra 
dim [su(N)\ = 2N(N - l)/2 + (N - 1) = (N 2 - 1) = d. 

In order to reverse engineer an exact quantum dynam- 
ics in su(N), one can follow the following five-step pre- 
scription: (i) First, choose r = N — 1 arbitrary frequency 
functions for the mutually commuting Cartan operators, 
u>j(t), with j = 1,2, ...,r. (ii) Second, choose N(N — 
l)/2 arbitrary two-dimensional vector fields, Xv{t), a n d 
use them to build an SU(N) rotation as a product of 
N(N - l)/2 non-commuting SU(2) rotations, U N (t) = 
Yl„ {"Tv-fi cos Xv{t) -i{n v (t) -T v )smxi,(t)}, c.f. Eq. (|37|. 



(iii) Determine the generators of the exact solution from 



<S> a (t)T N [J a ] = £ co 3 (t)W N (t)T N [hj]U N (t). Note that 

the latter is a doable straightforward calculation in the 
iV-dimensional representation described above, (iv) De- 
termine the unitary rotation and the diagonal opera- 



tor T>(t) — f ds exp < —is ^ LOj(t)Tjy2_ 1 [hj] > , which 

{ 3 =1 J 

are expressed in terms of the simple universal func- 



tion w(z) in Eq. (42) but involving some non-universal 
combinations of w's. (v) Finally evaluate b(t) = 
U]^2_ 1 (t)V{t)U N 2_ 1 (t)&{t), which therefore provides a 
precise driving field needed to perform the desired evo- 
lution. 
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Let us mention here that the su(4) case may be of 
interest to quantum computing. A particular realization 
of this 15-dimensional algebra is a pair of coupled half- 
integer spins with the most generic Hamiltonian of the 
form 

H sm (t) = bi(t) • &i + ba(t) ■ & 2 J afi {t)at,c*a%f), 

a,/3 

(44) 

where & a (i) are the standard Pauli matrices and we allow 
all the 15 parameters to be functions of time. The full 
Hilbert space of the problem is spanned by the wave- 
functions of the form, \tp) = c x |ft) + c 2 |tl) + c 3 |4-f) + c 4 l-U 
) . The normalized states reside on the seven-dimensional 
sphere and include, in particular, entangled states. Using 
the proposed five-step procedures, one may engineer time 
pulses that "untie" the entangled states back into pure 
states and vice versa thereby untying the corresponding 
second Hopf fibration, S 7 / S 3 — S 4 , which was proposed 
to describe two-spin entanglement [71] [72] ■ 



F. Two remarks on the global exponential 
representation 

1. Square-root and fractional powers of the evolution matrix 

If the dynamic symmetry group can be represented 
globally as an exponential of algebra elements, this rep- 
resentation leads to a very curious "natural" definition 
of matrix powers, as follows. Assume that G.4 is a dy- 
namic group and consider its element U = exp [— i*f> ■ J] . 

Define an arbitrary real (a € M) power of U as follows 



U a = exp [-ia* -J]eG A . 



(45) 



For example a square root of an SU (2) matrix (in two- 
dimensional representation), U = e5*' CT , would be sim- 
ply \fu = ei*°. 

It appears that using such fractional powers one 
can construct mixed Schrodingcr-Hciscnberg represen- 
tations of quantum evolution (which may conceivably 
be used for derivation of generalized "rotating-wave 
approximations"). Consider for example an average 
of operator products {not a Lie-algebraic invariant), 
C ab (t) = (V(t)\J a j b \y(t)) with |tf(t)) = U t \iPo) = 
exp — z<& • J \i/jq). We can write it as follows 



c ab {t) = (vj(o)\u t *u t ~ 2 Ja\/u t u t *J b \jUt\ju t \i> ) 

(46) 

or equivalently as 

c ab (t) = (^cv^^ija/^^jct/^^^d/^^^ (47) 



where |^ (1/2) (i)) = exp 



|^o> and ti 1/2 \t) = 



e2 * J J a ar e the wave-function and operator in such a 



While it is admittedly unclear whether the exotic 
mixed representations and our ability to raise group ele- 
ments to an arbitrary power are useful for practical calcu- 
lations, the existence of this operation brings up a more 
serious question about its multivaluedness and most im- 
portantly about the "analytical properties" of the global 
exponential map. This leads to the following section. 



2. Parameterization problem 

Here we point out an important problem of 
parametrization of the dynamic symmetry group that 
is present in both the DSBE "single-exponential ap- 
proach" [55] advocated here and product representation 
methods routinely used in the mathematical and quan- 
tum control literature [501 EH EE], he., where U{t) — 
d 



n« 

a=l 



,-i X a(t)J a 



(no summation). Namely, the problem is 



to determine a domain where the dual generators, «&(£) 
[or x{t)]i are actually defined. A naive view is that 
they can be taken arbitrary d-dimensional vector fields 
in a one-to-one correspondence with the driving field, 
b(t), which at any given time, t, is a vector of the d- 
dimensional vector space, M. d , associated with the un- 
derlying (i-dimensional algebra, A. This naive approach 
however can not be entirely correct, because the gen- 
erators are supposed to parameterize group elements, 



and the group G4 



is in most cases quite dif- 



ferent from and in some intuitive sense "smaller than" 
the Euclidean space, M. d ~ A*. For example, the su(2) 
algebra is just M 3 as a vector space, while its cover- 
ing group, SU(2), is topologically a three-dimensional 
sphere. This entire SU (2) group is reproduced in the two- 
dimensional defining representation, via the elementary 
identity: U[&] = e - **-* = I 2 cos [f ] -i(n-&) sin [f ] . 
Clearly, if $1 = $ 2 +47rfc and «i = n 2 or if $1 = Airk— 4> 2 
and rii = —n 2 (with k € Z + being positive integer), both 
$1 and <& 2 give rise to the same actual group element. 
Interestingly, these twisted conditions are much reminis- 
cent to the Hopf fibration, S 2 = S 3 /S 1 , with the modulus 
of dual generator, $ = |3>| to be associated with S* 1 and 
its "direction," n with S 2 . This example suggests that in 
fact proper parameterization of the dynamic groups may 
be far from trivial. 

Let us now recall that what seems to be a standard ap- 
proach in the mathematical literature to relating the Lie 
group and its Lie algebra is to define the elements of the 
latter as equivalence classes of curves passing through the 
identity in the former. While the author is not qualified 
to provide any valuable insight into the mathematical 
aspects of this old mathematical construction, we would 
like to speculate that in the context of quantum dynami- 
cal systems, a reverse procedure, if at all possible, would 
have been of great interest in that the algebra rather 
than the dynamic group is suggested as a primary struc- 
ture. An interesting mathematical question is whether 



mixed l/2-Heisenberg-l/2-Schrodinger representation. one can define group elements of Gyi 



J A 



as equivalence 
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classes of the algebra elements instead along the lines 
of the Hadamard construction in [73]. Let us conclude 
this somewhat obscure mathematical discussion, by not- 
ing that the formal question of parametrization of dual 
generators may in fact be of importance to physics, where 
proper parametrization of the dual generators may cir- 
cumvent the sign problem in certain interacting theories, 
see, e.g., Sec. | VII B | where a dual Hubbard-Strtonovich 
approach is presented for a quantum spin model. 



QUANTUM-TO-CLASSICAL 
CORRESPONDENCE 



We have been calling Eq. (151 "the dual Schrodinger- 



Bloch equation," but did not explain why we have cho- 
sen to associate Bloch's name with the equation that 
seemingly focuses on quantum dynamical systems only. 
We show below that the equation is in fact equally ap- 
plicable to describe the corresponding classical evolu- 
tion. Note that we have been using the units where 
the Planck constant is set to one, H = 1. Had we 
not chosen this convention for the units, we nevertheless 
would still have found that in most sensible physically- 
motivated cases, the DSBE does not actually include the 
Planck constant anyway and therefore is unchanged in 
the semiclassical limit, h — > 0. Also, it has been em- 
phasized that the Lie algebraic approach is by construc- 
tion representation-invariant, and hence all results apply 
equally well to large-dimensional representations, includ- 
ing the limit L — > oo, usually associated with the quasi- 
classical approximation. These arguments suggest, and 
the simple discussion below explicitly shows, that the 



DSBE equation (15 1 is indeed as applicable to describe 



classical dynamical systems as it is to describe the corre- 
sponding unitary quantum evolution. 



A. Generalized Bloch equations 

In a typical non-equilibrium quantum-mechanical 
problem with a well-defined Hilbert space, Hil(L), we 
are interested in calculating observables such as averages 



M(t) = (^|J(i)|Vo> 



(48) 



and correlation functions. E.g., a general two-point cor- 
relator can be presented as 

C{t u t 2 ) = ^c ab (< 1 ,t 2 )(^ |J a (i 1 )J b (< 2 )|Vo), (49) 

where |^o) € Hil(L) is a wave-function that describes 
an initial state of the actual physical system and J(i) = 
W (t)JU(t) are Heisenberg operators acting in Hil(L). 
Note that these Heisenberg operators belong to the alge- 
bra, because 



J a (t) = exp 



-iad 



*(t)-J 



which involves nothing but commutators. On the other 
hand, the observables such as the matrix elements (48) 



depend on the Hilbert space and the initial conditions 
and are not algebraic invariants. However, the amount 
of information we "need from the Hilbert space" is 



very small for the purpose of calculating (48 1 and re- 



duces to the initial conditions for the averages, M(0) = 
("Pol Jl^o)- I n fact, one can specify arbitrary initial con- 
ditions without reference to any Hilbert space (or quan- 
tum mechanics all together), with the caveat that an ar- 
bitrary classical choice may fail to satisfy quantization 
requirements, which automatically arise in the quantum- 
mechanical formulation. Given a set of initial conditions, 
M(0), the full time-dependence follows 



M(t) = i?[*(<)]M(0), 



(51) 



where R[&(t)} 



-F[#(t)] 



are the matrices uniquely de- 



termined by the solutions to the DSBE and therefore 
they are Lie-algebraic invariants. As to correlations func- 
tions, only averages of the commutators are expressed in 
terms of Lie- algebraic invariants in the same way as M(i) . 
E.g., a correlator involving a commutator of two opera- 

tors C aita2 (ti,t 2 ) = -i(ipo\ J ai (h), Ja 2 {h) \i>o) evolves 
into 

C ai , a2 (t u t 2 ) = R a b i i(t 1 )R£(t 2 )f bib c 2 M c (0), (52) 



where f bib2 are the structure constants (|4| of the under- 
lying algebra, A, and the tensor R(t) is a Lie-algebraic 
invariant defined via (51). Therefore to find the corre- 
lator, C aia2 , the only required piece of information that 
is not representation-invariant (and thus may depend on 
the Hilbert space) is the initial conditions, M(0). Let 
us note here that this not generally true for an arbitrary 
correlation function of type (49), which involves opera- 



tor products apart from "pure commutators," or more 
generally for any observable that is not a member of the 
algebra. This can be seen by assuming the contrary and 
finding a counter-example, e.g., in, perhaps, the simplest 
case of the two-dimensional representation of su(2) with 
the "correlator" (f |d^(t)|t) = 1, which is identically 
equal to one independently of evolution, while for all 
higher-dimensional representations, 7 1 | >2 [J Z ] is not pro- 
portional to the identity matrix and may exhibit a non- 
trivial time-dependence. 

Let us now confine ourselves to considering Lie- 
algebraic-invariant dynamics in the sense defined above, 
i.e., averages, correlators, Berry phase, etc. The averages 
can be described in terms of generalized Bloch equations 
that follow from the Schrodinger equation ([I]), or equiva- 
lently, from the familiar Heisenberg equations of motion 

i(t) =i[u{t), J(t)] , 

where the Hamiltonian, H(t) = b(i) • J and the structure 
constants @ of A yield 



J a = R a b (t) J b , (50) 



M b (t) = f ba c b a (t)M c (t), 



(53) 
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which have to be supplemented by the initial conditions, 



M(0). The generalized Bloch equations (53) are the ana- 



logue to the classical equations of motion for an arbi- 
trary Lie algebra. For some solvable algebras, such as 
the harmonic oscillator algebra, they reduce to the New- 
ton equations of motion in a Euclidean space (or more 
precisely to Hamilton's equations of motion). For the 
simple algebra, su(2), these are the usual Bloch equa- 
tion on the two-dimensional sphere, M = b(<) x M. In 
general, equations ( |53[ ) represent classical dynamics on a 
non-trivial manifold with constraints, .M, embedded into 
the d-dimensional Euclidean space. 

We have shown that the solutions to both quantum 
and classical dynamics are uniquely determined by the 
generators, & a (t). The classical solution is governed by 
the action of the operator, R, in Eqs. ( 51 ) or equivalently 
below 



M(f) 



= e -*<*> 



'M(O). 



(54) 



Recall that the classical averages form a d-dimensional 
vector that is acted upon by the (d x (i)-dimcnsional ma- 
trix, F [<&(£)] = <f>(t) ■ f from the adjoint representation 
of the underlying algebra. The Casimir invariants con- 
strain the vector of classical averages to lie within a clas- 
sical manifold, M.. In most cases of physical significance, 
this manifold is a symmetric Riemannian space with the 
metric that is uniquely determined by the Cartan-Killing 
product (16 1. Therefore, Eq. ([541 represents action of a 



symmetry group, G(M), on the classical manifold, M. 
Topological construction of the classical manifold as a 
coset representative of the dynamic symmetry group is 
discussed, e.g., in Ref. [18]. For completeness, we also re- 
view it in appendix [B] where it is argued that in the con- 
text of dynamical system at hand, the Perelomov's con- 
struction, which is slightly different from that of Gilmore 
in that it allows an arbitrary reference states to be used, 
is preferable and it is natural to associate the reference 
state with the initial condition, | V'o) 7 so that the classi- 
cal manifold is A4 = G a / H (\ipo)) , where H(\ipo)) is the 
maximum stability subgroup of the initial state. 

Note also, that the classical symmetry group is gener- 
ally expected to be smaller than the covering group gen- 
erated by the algebra, i.e., G{M) C G_4 = exp(i_4). E.g., 
in the previously considered case of the su(2) algebra [57], 
G(S 2 ) — SO(3) is the group of three-dimensional ro- 
tations to which SU(2) = G su (2) is a double-cover, 
SU(2)/SO(3)=Z 2 . 

Concluding this section, we reiterate its main conclu- 
sion that quantum-to-classical correspondence manifests 
itself in that the dynamics of classical averages ( 54 ) and 



correlation functions ( 52 ) are governed by the same gen- 



erators, 3?(i), that determine evolution of the quantum 
wave-function, see Eqs. ((8|) and (12 1, and satisfy the dual 



Schrodinger-Bloch equations ( 15 ) 



B. Example of quantum-to-classical 
correspondence for a linearly-driven 
Jaynes-Cummings model 

Here we provide an explicit simple example that illus- 
trates quantum-to-classical correspondence and demon- 
strates the Lie-algebraic approach in action. Consider 
the driven Jaynes-Cummings Hamiltonian: 



H lc (t) = A b (t)Pb+A s (t)'- 



2 



Pa- + ba A 



(55) 



Here P jb are canonical boson creation/annihilation op- 
erators, a z and a± = (a x ± ia y ) are the standard Pauli 
matrices, and A b (t), A s (i), and g(t) are arbitrary func- 
tions of time. As the title of this section indicates, 
we will be specifically interested in the linearly-driven 
Landau-Zener problem with A s (t) = —A b (t) = Xt and 
g(t) = g = const. 

First note that the operator N = Pb + ^f- commutes 
with the Hamiltonian and any other operator present in 
Eq. ([55]). I.e., 



N, Pb 



N, a z 



N. Pa- 



N, ba. 



0. 



(56) 



Therefore, we can treat the operator N as a c-number. 
Second, let us also define the following operators 



2V2F+I 



ba+, 



1 

2y/2N + 1 



Pa. 



and 



t a =Pb-N = 



These operators span the su(2) algebra, i.e. 



]+, £_ 



2S 3 and 



£a, Sh 



(57) 



(58) 



(59) 



(60) 



The Hamiltonian ( 55 ) becomes 
Hj C (t) = A b (t)N + A 3 (t)t 3 + g(t) (s + + S_ ) . ( (i 1 l 



with A 3 (t) = [A s (t) - Ab(t)] and g{t) = y/2N+lg(t). 
The first Abelian term in Eq. (61) gives rise to a triv- 
ial phase dynamics and can be factorized away. In the 
specific case of linear drive, we end up with the Landau- 
Zener Hamiltonian for %y, = 'Hjc(i) — Ab(t)N: 



= 2AtE 3 



9 



(62) 



The evolution operator ~ Since the algebra of E- 
operators is isomorphic to su(2), we can look for solu- 
tion in any faithful representation, including the two- 
dimensional spin- 1/2 representation. This maps the 
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original problem onto the exactly-solvable Landau-Zener 
problem in its canonic form |74) : 



= [\ta z +ga x ]V(t), 



(63) 



where ^f(t) is a two-component spinor. Note that the 
symmetry of Eq. ( 63 ) with respect to time-reversal de- 
mands that if tyx(tj~= ("') is a particular solution, then 

^2 (t) = (^.') is also a solution. If an arbitrary par- 
ticular solution is known, then a full S'-matrix can be 
constructed as follows EH 



in terms of (b(t)) and (<r(t)). The latter representation 
was used [76] in the context of the Dicke model (which 
in effect is a non-trivial generalization of the Jaynes- 
Cummings model for a higher-dimensional representation 
of spin). The "semi-classical approximation" used there 
involved replacing the average of operator products with 
a product of averages. For certain driven dynamics, this 
approximation can be shown to become reliable when the 
dimensionality of the representation (or equivalently the 
number of spins in the Dicke model) becomes large, but 
it is the least reliable in the opposite limit of the Jaynes- 
Cummings model. 



i'O 



_ / u u t - v v_ t u v_ t - v u t 

2 I U* V t - V U*^ t U U*_ t - V^Vt 



(64) 

where Vl/i(0) = ("°) is an initial condition that the known 
solution satisfies. 

A particular exact solution was derived, e.g., by Lim 
and Berry [75] as follows: 



*i(t) 



D u fo(t)] 



(65) 



where 7 = g 2 /A, v = ij/2, and D u (z) are known 
parabolic cylinder functions with well-studied asymp- 
totes. They give rise, in particular, to the standard 
Landau-Zener exponential transition rate as well as to 
non- analytic corrections to it, as discussed in the Ap- 
pendix of Ref. 175 

5J, and the identity, Ui^z{t) = 



From Eqs. (|64| 

3*Lz(t)'CT 



e 2 - = i 2 cos (— —i («LZ • o-) sin (— J^), we can 

find the representation-invariant and Hilbert-space inde- 
pendent dual generators. The evolution operator for the 
original Jaynes-Cummings problem therefore reads (here 
we restored the trivial phase dynamics): 



JO 



exp 



A 



t 2 N - 



#LZ(*) ' S 



(66) 



Note that short of exotic initial conditions with entangled 
spin-boson states, simple initial states (e.g., "spin- up" 
at t = 0) will evolve in exactly the same way as they 
would in a pure spin- 1/2 problem with the same driving 
field and the presence of the bosonic bath will manifest 
itself only in a rescaling of the Landau-Zener parameter. 
Another way to make the same statement is to focus 
on evolution of the classical averages, which occurs via 
action of the 5*0(3) rotation group 



[2N - 
i[2N 



i]-*<St*_ 

- 1]-§(S(T + 



-ba+) 
ft a-) 



= M{t) = e^ z(t) -fM(0), 
(67) 

where the (3 x 3) finite-rotation matrix, -F S u(2)[3?lz(£)], 
is given by Eq. ( [28] ) . Note that the classical theory is for- 
mulated here in terms of the averages of products of the 
bosonic and spin operators, as opposed to a formulation 



VI. LIE-ALGEBRAIC APPROACH TO 
INTERACTING MANY-PARTICLE QUANTUM 
MODELS 



Lie-algebraic formulation of a quantum lattice 
model 



So far we have focused on single-particle quantum dy- 
namical problems for a Hamiltonian linear in algebra 
generators. We now address a more general Hamilto- 
nian associated with an interacting many-body quantum 
system. The general system we have in mind is a lat- 
tice model with pair-wise interactions between particles 
associated with the sites. In the Lie-algebraic sense, 
it is equivalent to having a representation of a finite- 
dimensional algebra associated with each site (e.g., Fu 
for bosons, u(2) for spins, u(3) for "quarks," etc.). 

Define a discrete set, S, e.g. corresponding to a lat- 
tice in physical space. Assign to each site i € S, a finite 
dimensional Lie algebra, A(i) = span { Jo(i), ■ ■ ■ , ^(i)} 
and its L^-dimcnsional representation. Define an inter- 
acting Hamiltonian as a form bilinear in generators of the 
single-particle algebras. In what follows, we shall focus 
on the case where the single-particle entities on all sites 
are identical to each other, i.e., Vi E S, A(i) ~ A and 
the corresponding representations are the same L(i) = L. 
This is not the most general case that can be considered 
(e.g., the text-book Dicke model discussed in the previ- 
ous section does not belong to this category as it involves 
two different algebras, f} 4 and u(2), on two "sites"). Let 
us also require that the algebra A includes an Abelian 
u(l) component, Jo (i.e., [Jo, Ja] = 0), which always 
can be incorporated into the theory. In what follows, we 
identify it with the unity operator, J a (i)Jo(j) = Ja{i)- 
Note that for the interacting problem, it is important 
to specify the physical dimension of the single-particle 
Hilbert space as the many-particle algebra generated by 
the interacting Hamiltonian will generally include pow- 
ers, J", which acquire a precise meaning only given a 
representation. With this, let us write down the inter- 
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acting Hamiltonian in the following form 



7i = £5> a (W) 

ies 0=1 



1 

2^ 



J2 Y.^o)Ui)Uj), 

a— 1 

where fis is the number of sites in the discrete set, S, 
H a are real parameters that may correspond to external 
fields or a chemical potential, and X a (i,j) = X a (j,i) is a 
symmetric interaction matrix. Note that in Eq. (68 1, 
we have excluded from consideration interactions be- 
tween different elements in the algebra, but even with 
this simplifying assumption the model ( 68 ) still remains 



very general and includes a huge number of condensed 
matter Hamiltonians (excluding certain topological lat- 
tice models). For example, the hopping terms of the 



Bose-Hubbard model can be written as, b[b 2 + b 2 bi = 
X1X2+P1P2, where £ = (b + w)/yf2 and p = i(tf — b)/y/2 
on any site. Note that the usual density-density inter- 
actions, b\b\b\b2, also are bi-linear in terms of genera- 
tors b^b from [54. In this Lie-algebraic sense, there are 
not considered any more complicated than the hopping 
terms, eventhough writing them down involves an ex- 
pression quartic in terms of boson creation/annihilation 
operators. 

There are two ways one can proceed along in con- 
sidering the Hamiltonian (68). A first approach is to 
define new composite operators XfJ' = J a (i)Jb(i), and 
take all non-zero terms in the Hamiltonian or their spe- 
cial linear combinations and commute them with one 
another, (^ad^X'^j , continuing doing so until the al- 
gebra closes. For any finite-dimensional representation 
of single-particle algebra, A, it is certainly guaranteed 
to happen because any matrix satisfies its characteris- 
tic equation and so the number of linear-independent 
powers of X n is limited (note however that for infinite- 
dimensional representations, e.g., for the boson case with 
no restriction on their number on a site, it is not gener- 
ally so). In any case, the many-particle Hamiltonian (68) 



then can be viewed as a linear form in and a member 
of the new, possibly huge, many-particle algebra, A- m t, 
that gets generated this way. Such an approach may be 
a sensible thing to pursue if the many-particle algebra 
formulated in terms of some physically-relevant genera- 
tors closes into a tractable form or factorizes into sim- 
pler components. This would not be unexpected in inte- 
grable cases, or when the rank of the many-particle alge- 
bra (or equivalently the number of conserved quantities 
or Casimir operators) grows steadily with the "system 
size." According to the preceding "single-particle" sec- 
tion, as long as the many-particle algebra remains finite- 
dimensional, the classical and quantum dynamics (formu- 
lated in terms of the new generators of A ln t) are related to 
each other, and therefore classical integrability and quan- 
tum integrability are in a one-to-one correspondence as 
well, in the sense that the dynamics of the averages in 

the many-particle algebra, (X(t)\, determines the dy- 



namics of the operators and their commutators. Again, 
the existence of such correspondence may be established 



via the dual Schrodinger-Bloch equations (15) for A 



although its practical application may be challenging 



B. Hubbard-Stratonovich Dynamical Systems 

However, a general case, relevant to most "real-life 
Hamiltonians," is that the many-particle algebra "ex- 
plodes" fast and often becomes infinite-dimensional. Also 
the integrability status of such generic models is rarely 
known. Having this sort of a case in mind, we proceed via 
a different route, generalizing a construction of Ref. [54] . 
We concentrate on the thermodynamic problem first, i.e., 
on a time-independent interacting Hamiltonian ( 68 ) . Wc 



assume the applicability of a path-integral representation 
for the following quantity 



Tre 



-Hu 



(3 + it, 



(69) 



where the trace is over the entire L^ s -dimensional 
Hilbert space of the many-particle problem. Notice that 
we allow the parameter u in Eq. ( 69 ) to be a complex 
number such that u = j3 = 1 [T corresponds to the parti- 
tion function, while Z(u — it) reflects the property of a 
unitary evolution |77] (e.g., describing a free evolution af- 
ter a quench in system parameters) . Now assume that the 
usual path-integral representation of the partition func- 
tion, Z, applies, i.e., 



Z(u) 



M fc (0)=M fc (u) 



Y[V[M k ( T )}e 

fces 



\ rx.p^ - I dTH[M k (T)}j , (70) 

where Sherry [M^] is a topological term on a site k, en- 
forcing proper commutation relations for single-particle 
generators, and H [Mfe(r)] is a Hamiltonian (68) with all 
operators replaced with the corresponding coherent-state 
vector fields. The integration over r in Eq. ( 70 ) is along 
a straight line in the complex plane, t = r' + ir" , that 
goes from to u = j3 + it and the coherent-state trajec- 
tories are required to form loops, M(0) = M(u). Note 
that the existence and validity of such a path-integral 
representation is a "big if" from a formal mathematical 
standpoint. One problem related to the previous dis- 
cussion is the choice of a classical manifold, M., where 
the loop-trajectories reside. The single-particle construc- 
tion does not necessarily apply here, because it relied on 
single-particle pure states, while here entangled states are 
possible and generally unavoidable. This is however not 
a concern if group action covers the entire Hilbert space 
(e.g., spin-1/2 systems) and also some version of the path 
integral construction is more likely to remain reliable 
in cases where the classical manifold is representation- 
independent (e.g., the lattice bosonic models associated 
with the solvable single-particle algebra, (54). At this 
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point, we just assume that this is so and that classical 
trajectories remain in a singly-connected classical man- 
ifold M(t) C A4. Now let us focus on the interaction 
term in the action 



exp 



p 



We can decouple the interaction terms for each non-zero 
^ a (hj) 7^ 0, using the following identity 



-z u XM 1 M 2 _ 



dA _a? 

-e 2 i A 



+ ii„A(M 1 2 +Mf ) 



exp { V-^sgnAA (M x + M 2 ) } , (72) 



where A € K and z u — u/\u\ = {(i+it)/ \J j3 2 + t 2 is a unit 
vector in the complex plane associated with the choice of 
the time-integration contour in Eqs. (69 1 and (70 1 (e.g., 
for thermodynamics, it is just equal to one, zp = 1). If we 



implement this transform ( 72 1 in the path-integral ( 71 1 , 



it renormalizes the on-site interactions as follows, 



\ a (i,i)^\ a (i,i) = A a (i,i)- \ A a ( 



(73) 



[where the summation is over all interacting links con- 
taining a given site, but it does not generate any 
other terms. We now can introduce another set of 
Hubbard-Stratonovich fields on sites, Bi(r), to decouple 
the remaining on-site interactions to arrive to a repre- 
sentation in which for each particular realization of the 
Hubbard-Stratonovich fields, Ai(t) and i?i(r), the re- 
maining path integrals over the coherent states decou- 
ple on different sites and correspond to a single-particle 
driven dynamics, where the time-evolution in the algebra 
is governed by the corresponding Hubbard-Stratonovich 
fields. Therefore, we can "roll back" the path-integral 
construction into the Hamiltonian formulation to find the 
following main general result 

oo 

Z(u)= J lipAfir)} II pBt(T)]Y[ Zi [A,B] 



a;/£liiiks 

x exp < — — / ar 



E 



a;«Gsitcs 
2 



\A?(t)\ 



E 



\B?(t) 
|A?I 



(74) 



where I denotes all lattice links I — (i,j) for which the 
interaction is non-zero, i is a site index, and 

Zi = Trpi(u) 

is the trace of the time-evolution matrix that governs 
the behavior of the single-particle quantum-dynamical 
system, where the driving field depends on a particular 



Hubbard-Stratonovich realization 

-d T p(T)p- 1 (r) = H l (T) 

p(o) = i 



(75) 



We will call this system of differential equations, the 
Hubbard-Stratonovich dynamical system. The corre- 
sponding Hamiltonian is generally non-Hermitian and 
has the form 



Hi{r) 



d 

E 

a=l 



p a (t)-Ft(r)-J-z u XtBUr)\J a 



with F?(t) 



(76) 

^2 v~5uAf^f (t) and A defined in 

l=(i,j) 

Eq. ( 73 1 . Note that for a generic interacting Lie-algebraic 
model of type (68), the time-evolution, p, is not neces- 



sarily unitary, nor it is associated with any real form. 
From the Lie-algebraic point of view, the main difference 
between the Schrodinger equation Eq. Q and Eq. (75) 
is that the latter is generally to be implemented in a 
complex extension of the underlying single-particle alge- 
bra. This leads to a number of catastrophic complica- 
tions such as a break-down of the global exponentiation 
conjecture, which relied on the compactness of the under- 
lying dynamic group (however, as noted in Appendix |A 4| 
certain non-compact dynamical groups still allow global 
exponential map, see also next section VII A[ ). In the gen- 
eral case, we are forced to look for solution in a product 
form, see, e.g., Wei and Norman [601 161 j . or in a mixed 



form of type, p(r) = e 



-R(t 



i#(t)-j^ w hich may be 



"admissible" in certain cases (e.g., any SL(2,C) matrix 
can be represented in this form). 



VII. EXAMPLES OF THE LIE-ALGEBRAIC 
APPROACH TO MANY-PARTICLE LATTICE 
MODELS 

A. Non-unitary quantum evolution on the example 
of the Bose-Hubbard model 



Consider the Bose-Hubbard model [781480] on a two- 
dimensional simple square lattice in equilibrium, 



BH 



E 



Mo 



hi + u fij (hj - 1) , 



where t is the hopping amplitude, po is the bare chemical 
potential, U is the on-site repulsion, the indices i and j 
label sites of the square lattice, and I = (i,j) label its 
links. The bosonic creation/annihilation operators sat- 



b 



S%jl, 



isfy the familiar commutation relations, 

hi = b\bi is the density operator. On every site, r^, b\ 
b t , and 1 span the harmonic oscillator algebra. 
The partition function of the model is 



z bh = / n [^m] e 



-S[b,b] 



(78) 
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where the action is 



Sbh = ~J dr\^2 [^( T ) ( d r ~ M) h{r) + U \bi(r)\ 
o 

-tJ^HrMr)}. (79) 

Here we used the prescription for operator-ordering of 
Ref. [81] . which was suggested in order to reproduce the 
correct result in the single-site toy model. Therefore, 
p = p + U in Eq. (1791). 



The conventional approach is to use Bogoliubov mean- 
field or Gross-Pitaevskii analysis that lead to a reliable 
mean-field description of the superfluid phase and classi- 
cal fluctuations deep in the phase. This description how- 
ever becomes unreliable on the Mott insulating side and 
in the vicinity of the supcrfluid-to-Mott-insulator transi- 
tion as well as in the classical high-temperature liquid. 
To build an unbiased approach, we notice that from the 
Lie-algebraic point of view, the interaction terms are not 
any more complicated than the hopping terms, as both 
these terms represent bilinear forms in terms of elements 
of the harmonic operator algebra. In accordance with 
the general "prescription" from the previous section, we 
decouple the hopping terms first using the following al- 
gebraic identity 

£ T t(6 i b J +6 J fc.) _ f d 2 ai -£ T [||q i | 2 + (aiai + B i a,)-t("i+« ; )] 



e T n/t 

where I = denotes a nearest-neighbor link (the links 
themselves form a square lattice), e T corresponds to a 
single imaginary time-step in the truncated path integral, 
and 



B,(t) = & 4 (t)-& j (t), 



(80) 



where we adopt the convention that for the horizontal 
links, I — + x), the plus sign in Eq. (80) corresponds 
to the i-th site on the left, and for the vertical links, 
I = (i,i + y), the plus sign corresponds to the z-th site in 
the bottom. This yields 



-<BH 



n [©\(t)] n [^m] 



ri exp 




h{r) (d T - p) bi(r) 



+ {^(tM^ + c. c.} +U\b l (T)\ 4 \ J, (81) 
where p — po + U — 4t and 

Fi{r) = curLa(r) = a^ ti+it ) - a(,_£,j) +a^ >i+ ^) - 

(82) 

is a curl of the Hubbard-Stratonovich field on the site i, 
which is a plaquette of the link lattice. 



1. Hubbard-Stratonovich dynamical system for the 
harmonic oscillator algebra, f)4 



Let us now follow the general prescription of Sec. |VI B| 
and decouple the on-site interaction as well with another 
set of real Hubbard-Stratonovich-fields to arrive to the 
following general expression 



-BH 



n^M]n^M]^ i? ' drMr)|2 



(83) 



where the imaginary-time evolution operator, Pi(r), sat- 
isfies the following non-Hermitian Schrodinger equation, 

-d T pi(r) = {{iV l (T)-f t ]tfb+ [^(r)S + 7i(r)6t]}^(r), 

(84) 

with the initial condition, Pi(0) — 1. Note that there 
is known to exist an ambiguity in selecting a channel in 
the Hubbard-Stratonovich transformation as discussed in 
Refs. [82] , [83] , and [84] . I our case this subtlety manifests 
itself in that we can decompose the on-site term into 
linear combinations of either b'w and bb or b^b, which 
gives way to different single-particle algebras. 

We choose the latter decoupling, which leads to the 
action in Eq. (83) that is not only local, but also corre- 
sponds to the (non-unitary) dynamics in the f)4 solvable 
algebra. Using the method introduced in Sec. |IV A[ we 
find the exact solution to Eq. (84) as follows 



Pi(r) = exp -Wi(r)Pb 



(85) 



where Wi (r) 



-pr + i J dsVi(s), r)+,i( T ) 



fd S F t ( S )e w ^ s \ n- ti {T) =JdsFi(s)e- w *W, and 7l (r) 
o o 

fdsj d« / Ji(a)Ji(fl / )ex P [W , i (fl') - Wi{s)}. 





To analyze the path integral ( 83 ) , we need to calculate 



the "dynamic partition function,' 



n=0 



(86) 

/ a \ 71 

where \n) = (6m /Vn\\0) arc the usual harmonic os- 
cillator eigenstates. As long as this series is convergent 
(the interaction U is not too large), this trace can be 
evaluated exactly utilizing the Glauber coherent states 
and the result reads 



aim 



exp 



1 - e-WiiP) 



(87) 
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The exact non-perturbative action is therefore given 
by the sum of these terms, 



S C B = 



ies 



In 



where the functions involved in Eq. ( 88 1 arc defined below 



Eq. (85) and we recall that ^(t) = curler), with ai(r) 
being a complex field defined on the links, V^(t) is a real 
field defined on the sites. The exact ■partition function is 
a path integral over both ai (r) and V, (r) , 



-"BH 



Yl^a^nmwe- 



[a,V] 



fl 



(89) 



We emphasize that the full action (88) contains simple 



Riemann integrals over imaginary time. Also, notice that 
had the field, V$(t), been absent, the partition function 
(871 would have been manifestly real and positive def- 



inite and therefore it would not exhibit any sign prob- 
lems. The presence of this field complicates matters in 
the quantum limit. However, since the main part of the 
action, S e g[a, V], has been evaluated into an explicit form 
and the other two terms are comparatively trivial, the nu- 
merical analysis of the action ( 89 1 may be an alternative 



to the conventional Monte-Carlo technique. 

Also note that the fields ai and Vi can be shown to 
be related to the superfluid and Mott order parameter 
correspondingly. Therefore, the exact non-linear action 
( [83] ) should contain a theory of the two competing or- 
ders, as well as a complete description of the superfluid- 
to-Mott-insulator quantum phase transition. However, if 
the interaction becomes relatively strong, the series ( 86 ) 
becomes manifestly divergent |85j and its careful resum- 
mation is required to get a meaningful result [86] . These 
complications are not unexpected because we know that 
the U — > oo "hardcore boson" limit gives rise to the 
XY quantum spin model |87j formulated in the com- 
plex extension of su(2), and the corresponding sl(2,C) 
Hubbard-Stratonovich dynamical system is not exactly 
solvable for an arbitrary quantum fluctuation field |54j . 
Note here that interestingly, the hardcore-boson limit can 
be viewed mathematically as a contraction of the infinite- 
dimensional single-site Bose-Hubbard algebra, 



span <n,n 



(90) 



into su(2) as U — > oo. Clearly, the dual description of 
the corresponding ATY-model should not include fluctu- 
ations of the site-field, Vi(r), which will implement the 
no-double-occupancy constraint and reduce Eq. (81) to 



the dual partition function of the following type, (here 
by "dual," we imply an unusual form of duality in lattice 
link-space, c.f., Ref. [88 , which is also different from dual- 



ity in the preceding sections and Sec. VII B which refers 



to the imaginary-time domain or equivalently quantum 
fluctuations) : 



■>XY 



\[ [V 2 ai{T)} e 'o 



(91) 



[J 1 + e" M0/3 exp 



p 

2Re J dr curl; a(r) 



Here, /iq is related to a uniform Zeeman magnetic field 
in the z-direction in the spin formulation and the curl of 
the link fields, a ; (r) is defined in Eq. (82). 



2. Hubbard-Stratonovich dynamical system for the 
infinite- dimensional single-site Bose-Hubbard algebra, f) 2 

Finally, let us present yet another approach to ana- 
lyze Bose-Hubbard model, which appears to be especially 
useful for constructing various doping-dependent mean- 
field approximations of the theory [86]. Notice, that ac- 
tion (81) was already local in site- fields and we intro- 



duced the additional on-site Hubbard-Stratonovich fields 
in order to reduce the problem to a tractable dynamical 
system for the solvable algebra, t)4. Alternatively, we can 
bypass this step and present the partition function in the 
following form instead 



•-BH 



n[v\(r)} e ^?^ r|ai(r)|2 nTrA(/3){92) 



where now the "density matrix" is governed by the fol- 
lowing Hubbard-Stratonovich dynamical system 

- d T Vi(T) = [Utit - fail, + [>,(>)&« + .Fi(V)&t] } V^t). 

(93) 

Notice that the operators, |n|, hi, Sj, b\, i|, and their 

commutators do not close into a finite-dimensional al- 
gebra, but form an infinite-dimensional algebra instead. 
Eq. (93) does not appear to be solvable in the opera- 



tor form for an arbitrary time-dependence of its coeffi- 
cients. However, in the regime of strong interactions and 
on Mott insulating side, we can approximate the prob- 
lematic interaction term as follows, nf « hi(rii(r)) , which 
leads to a tractable self-consistent theory that retains 
quantum fluctuations. An analysis of the Bose-Hubbard 
model in the non-perturbative regime, including a care- 
ful analysis of its hardcore spin limit, will be presented 
elsewhere 186 . 



B. Dual action for a maximally- frustrated quantum 
spin model 

The single-particle part of this paper has focused on 
the dual approach to real-time dynamics for closed quan- 
tum systems in which the trajectory in the algebra, i.e., 
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the driving field, b(t), was related to the trajectory in the 
compact dynamic group, *&(t). It is interesting to see if 
a similar dual approach can find its use in the many- 
particle context. The motivation for looking into this 
question is that a dual description gives the solution to 
quantum dynamics right away, i.e. it immediately deter- 
mines "the effective action," S e ff = — ^ In [Trp.j (/?)], and 

i 

circumvents the need to solve the corresponding differen- 
tial equations. This approach is especially attractive in 
the cases where we integrate over all possible trajectories 
in the algebra and hence we may expect that such inte- 
gration can be replaced by a path integral over the dual 
fields with readily available quantum dynamics. 

As noted in the end of Sec. |VIB| the direct map- 
ping of the dual quantum-mechanical description based 
on the global exponentiation conjecture is not always 
possible for a generic Hubbard-Stratonovich dynamical 
system that arises in the many-particle context. This 



is because the corresponding dynamical system ( 75 1 , 
id T p(T x = /3(t) • J, involves complex driving fields that 
double the dimensionality of the real algebra into its com- 
plex extension by adding non-compact generators. Due 
to the non-compactness of the corresponding dynamical 
group, there is a fundamental restriction that in this case 
we can not write down the non-unitary quantum evolu- 
tion operator as a single exponential of algebra gener- 
ators, i.e., generally, p(r) 7^ exp j — [R(t) + £<&(r)] ■ J j 
and we only know that a factorized version of the solution 
is available [501 EI] ■ 

In some interacting many-particle models however the 
unitarity of quantum dynamics persists for Hubbard- 
Stratonovich dynamical systems as well. A rule of thumb 
for the Hubbard-Stratonovich dynamical group to remain 
compact is that the interacting model should include only 
repulsion (in the particle language) or antiferromagnetic 
terms (in the spin language). To provide an example, 
consider a quantum spin model with Heisenberg interac- 
tions 



H 




(94) 



As evident from (94), the model assumes all spins inter 



acting with each other, and therefore the geometry and 
dimensionality of the lattice play no role. Note that if 
the number of spins is three, N = 3, the model (94) be- 



comes a canonic toy model to illustrate the phenomenon 
of frustration. For a larger number of spins, it remains 
frustrated to a maximum due to the infinite-range anti- 
ferromagnetic interactions. Let us stress that the goal of 
this section is not to solve the spin model ( 94 ), but rather 



to illustrate the key features of the dual approach, which 
manifest themselves in more complicated models as well. 
The thermodynamic partition function of the model 



( 94 ) reads in the path-integral representation: 



Z= n\[VU 3 {T)]e * ' , (95) 



where Af^(r) is the unit vector-field on the Bloch sphere 
for the i-th spin, f2i(r) is the corresponding solid angle, 
and 5*wz [M] is the Wess-Zumino topological term in the 
action, whose explicit form is well-known |10j , but we do 
not need it (note that our general construction relies only 
on the existence of a path-integral). Let us now decouple 
the interaction term via a single Hubbard-Stratonovich 
field, A(t), as follows (we omit below an insignificant 
overall constant): 



Z= jv'A^Ylpn^r)} 



- ■£ Swz[A^j]+i / drA(r)-M(r) 

xe 3 



(96) 



Again, we can "roll back" the partition function into its 
Hamiltonian formulation to arrive to 



Z = Jv 3 A( T )e ""^ (T) /[A(r)], (97) 

where N is the number of spins and 2 [A] = Tr/5(/3), 
where the matrix p(r) follows from the Hubbard- 
Stratonovich dynamical system 



id T p{r)p-\r) = §A(t) 
P(0) = h 



(98) 



Due to the imaginary constant in Eq. ( 98 1, the dynamical 



group for this system is SU(2) and therefore the prob- 
lem maps onto unitary spin dynamics. Using the results 
of Sec. |IVB[ we can write the solution in the dual form 
p(r) = exp [— k«&(t) • <t] , where the dual field satisfies 
the dual Schrodinger-Bloch equations pO) and thereby 
connects it to the Hubbard-Stratonovich field. Further- 
more, using the elementary properties of the Pauli ma- 
trices, we find the "dynamic partition function," z[A], in 
Eq. ((97} as follows: 



z[A] = 2 cos [$(/?) /2] 



(99) 



which takes on a very simple "boundary form" in the dual 
language. On the other hand, the Gaussian A 2 -term in 
Eq. (97) becomes non-trivial. 



If we now replace the integration over the original 
Hubbard-Stratonovich field with that over its dual coun- 
terpart, we can write symbolically the following "dual" 
expression for the partition function: 



(100) 



r / 

: exp dr [n 2 (T) sin 2 



*(r) 



where *(r) = $(r)n(r) with \n(r)\ = 1 and det (ff ) 
denotes a Jacobian of the transform from the original 
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Hubbard-Stratonovich field covering all possible trajec- 
tories in the Euclidean space, M 3 to the dual field cov- 
ering all trajectories in the group, SU(2) ~ S 3 . Let 
us reiterate that in doing this transformation, we "trade 
off" the simplicity of the Gaussian term in the direct 
approach of Eq. ( |97[ ) for the simplicity of the quantum 
evolution term, Tr p(fi). The price we pay for this how- 
ever is that the former, initially trivial Gaussian term 
that penalizes any fluctuations takes on the form of a 
non- linear functional, see the last term in Eq. (100 1. Re- 



markably this factor in new variables describes a free par- 
ticle moving in a three-dimensional sphere. Indeed, if we 
set n = (sin cos 0, sin(9sin</>, cos 6*) and w\ — rcos9, 
W2 = r sin 9 cos cj>, W3 = r sin # sin </> cos ($/2), and W4 = 
r sin sin sin ($/2), then the Lagrangian in Eq, 



becomes L = 



n 



: (r) sin 2 



2 



2 



(100) 



Y^i—i provided that r — 2/«JJ is held constant. 

Finally, let us mention a curious property of the dual 
representation (100): if the number of spins is even, all 



terms in the partition function ( 100 1 are manifestly real 
and positive definite. This is to be contrasted with the 



direct Hubbard-Stratonovich representation (97), which 



exhibits a "sign problem" because different trajectories of 
the field, A(r), give rise to fast-oscillating terms in the 
path-integral. This observation suggests that the dual 
approach proposed here may be a promising avenue to 
cure the sign problem in certain theories. 



VIII. SUMMARY 

In this paper, we have put together a general Lie- 
algebraic approach to analyze quantum dynamics. The 
punchline of the first part of the paper is that for a large 
class of non-equilibrium quantum systems, the notion of 
a Hilbert space, central to the conventional Schrodinger 
formulation, is somewhat of a red herring, as the choice 
of a specific representation hides a more primitive Lie- 
algebraic structure of the Hamiltonian. We suggest and 
demonstrate on a number of explicit examples that the 
dual Schrodinger-Bloch equation for the Hilbert-space- 
invariant generators is a conceptually viable and practi- 
cally useful alternative to the conventional approach. Let 
us emphasize that this and some other results and state- 
ments of Sees. |TT] and |III| are not new and they appear 
to be scattered over the literature [S5HS3 on quantum 
control theory [51], generalized coherent states [18], as 
well as in the mathematical literature on the Lie the- 
ory and dynamical systems [68 . In fact, the equation 
identical to what we called here "dual Schrodinger-Bloch 
equation" appears in the introductory chapters of many 
mathematical textbooks on Lie algebras [55] m the con- 
text of the derivation of the Baker-Campbcll-Hausdorff 
formula. However, this particular form ( 15 ) of the equa- 



tion is not normally assigned major importance, perhaps 
because the mathematical researches often tend to study 
Lie groups in their utmost generality and it has been long 



known since the late nineteenth century that for generic 
non-compact groups, a global exponential map gener- 
ally does not exist. Due to this well-known fact, there 
has been no reason to study the equations of motion for 
"global generators" in the case of a generic dynamical 
system. On the other hand, as emphasized in this paper, 
a large number of quantum-dynamical systems seem to 
give rise to dynamic groups with surjective exponential 
map, where quantum evolution operator allows global ex- 
ponentiation from the algebra. This observation seems to 
enhance the importance of the dual approach based on 
the Magnus representation [58| and in fact one may ar- 
gue that whether to take the Schrodinger equation in a 
Hilbert space in its canonical form ([I]) , or its dual Hilbert- 
space-invariant version ( 15 ) as a starting point becomes 



a question of convention and convenience in this case. 

We further argued in this paper, that the dual gen- 
erators, &(t), provide a transparent and simple way 
to establish a correspondence between unitary quan- 
tum evolution and deterministic dynamics of the cor- 
responding classical system. A concise summary of 
quantum-to-classical correspondence is as follows: Quan- 
tum dynamics maps an initial normalized state into a fi- 
nal normalized state in an L-dimcnsional Hilbert space, 



!*(*)} 



exp 



-i#(t) • J L |tf (0)) c 5 2L - X c Hil(L), 



via action of a dynamic group in the L-dimensional 
unitary representation (i.e., Jl above are L x L Her- 
mitian matrices). Classical dynamics maps the ini- 
tial coordinate-vector onto a final coordinate-vector in 
such a way that both lie within a classical manifold, 
A4, embedded in a c£-dimensional classical phase space, 



M(i) 



exp 



*(t) • f d M (0) C M C R d , and this 



dynamics occurs via action of a symmetry group in its 
d-dimensional adjoint representation (i.e., fd are d x d 
matrices with real entries). Quantum-to-classical corre- 
spondence is in that the group action in both cases is 
governed by the same dual generators, 3>(t), and thereby 
determines the solution to both the Schrodinger equation 
and the corresponding generalized Bloch equations. The 
simplicity of this formulation raises the question of how 
general it is and also concerns in relation to chaotic sys- 
tems, where quantum-to-classical correspondence is far 
from clear-cut. Note that our approach hinges on the ex- 
istence of a global exponential map, which is guaranteed 
to exist for a certain (limited) class of finite-dimensional 
algebras. The general conditions to have such a map ap- 
pears to remain an open mathematical question under 
debate and even less is known about infinite-dimensional 
algebras. Therefore, our construction does not directly 
apply to quantum models with infinite-dimensional alge- 
bras, nor to chaotic systems. 

Let us also point out that while the dual approach does 
not perform the miracle of solving differential equations 
that are otherwise unsolvable in the direct formulation, 
it does appear amazingly useful in constructing exactly- 
solvable non-linear models. Ref. [57] and Sec. IV B have 



demonstrated that the method can be used to reverse- 
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engineer a variety of exactly-solvable two-level-system 



dynamics and Sec. IV E put forward a constructive five- 
step procedure of reverse-engineering exact solutions in 
more complicated higher-dimensional groups. In fact, the 
inspection of the latter five-step procedure indicates that 
it should be possible to start with a rather generic poly- 
nomial P{u,t) = Y[ [oj — Wfc(i)] l/fc , which with some mi- 

k 

nor restrictions on the functions Wfc(i) and £ Z + , can 
be viewed as a secular equation of an algebra, and then 
build upon it a quantum dynamical system, which would 
be exactly-solvable by construction. 

While quantum-to-classical correspondence and the 
general structure of the theory in the single-particle case 
(finite-dimensional algebra) appear to be understood, it 
is far from being so in the case of interacting many- 
particle models. In Sec. | VI A and Ref. [54], we attempted 
to extend the Lie-algebraic view to a certain wide class 
of quantum lattice models (excluding various topolog- 
ical models, where interactions between different alge- 
bra generators, or equivalently generalized "spin-orbit- 
couplings," would be allowed). The key new step pro- 
posed here is to use a generalized Hubbard-Stratonovich 
decomposition [54 , which maps a large, often infinite- 
dimensional, interacting algebra into an ensemble of 
finite-dimensional single-algebra (think single-particle) 
Hubbard-Stratonovich dynamical systems. One major 
difference between our approach and that of conventional 
Hubbard-Stratonovich decomposition is that we focus on 
algebra generators instead of counting the number of 
creation/annihilation operators in the conventional Fock 
representation. As demonstrated in Sec. VII A this ap- 



proach is quite useful especially for lattice models based 
on solvable algebras (such as the Bose-Hubbard model 
built of harmonic oscillators on lattice sites), where an 
arbitrary dynamical system, unitary or non-unitary, is 
explicitly exactly-solvable. 

In conclusion, let us note that an important wide-open 
mathematical question in addressing the many-particle 
lattice models of type ( 68 1 is the applicability of the 



path-integral approach itself [90 . In particular, a recent 
Letter of Wilson and the author [5T] has demonstrated 
that the standard normal-ordered path-integral seems to 
break down even for a model as simple as the single-site 
Bose-Hubbard model and while a way out was found to 
correct the theory into a sensible form, the latter non- 
normal-ordered construction is not based on any solid 
mathematical arguments, and therefore a careful investi- 
gation of this controversy is clearly called for. 
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Appendix A: General Properties of Quantum 
Dynamics 

The explicit form of the Hilbert-space-invariant DSBE 
( [T5| ) depends critically on the structure constants of the 
underlying Lie algebra. In this appendix, we make a few 
general remarks about unitary evolution and the struc- 
ture of Eq. (15 1, which can be classified in a one-to-one 



correspondence with classification of Lie algebras. 



1. Decomposition of the algebra into a standard 
form 

There exist powerful and effective methods of decom- 
posing an arbitrary Lie algebra into a standard form [68] . 
The first step involves factorization of A into a solvable 
subalgebra, _4o, and a scmisimple component, -4 sem sim- 
The Levi decomposition |91j below exists for any finite- 
dimensional algebra with no exceptions 



A = An + A 



scm sim 



A + J2 A ^ 



(Al) 



SUCh that [^Ascmsimi ^semsim] C *A scms i m , [-^0; -^o] ^- Ao, 

and [-4 scms im, Aq] C Aq- By definition, the scmisim- 
ple component may be decomposed into a combination 
of mutually-commuting simple sub-albegras: -4 scm sim = 
^2v=i -A-v, where we N B is the number of such simple sub- 
algebras. Note that our ability to actually decompose 
the algebra into this form and also to extract the simple 
components is tied closely to our ability to solve a generic 
eigenvalue problem on the algebra, which corresponds in 
our case to finding the spectrum of a generic stationary 
Schrodinger equation. 



2. Factorization of the evolution operator for 
solvable algebras 

We recall that a solvable Lie algebra is spanned by 



the generators A$ — span 



(0) 



yx 6 Aq, ad^J^ ' = LaJb'' (nilpotcnt algebras in- 

b>a 

elude a strict inequality in the sum, i.e., b > a). In the 
language of the quantum dynamical problem defined by 
( II I, including its DSBE form ( 15 ), it implies the follow- 



{./! 

-b f(0) 



j(0)\ 



such that 



ing: The non-equilibrium dynamics in the solvable alge- 
bra reduces to a chain of ordinary differential equations 
that can be solved one-by-one using the graded struc- 
ture. A solution can be obtained via factorization of the 
evolution operator as follows. Consider the Hamiltonian, 
which is a trajectory in a solvable algebra, Aq, i-Lo{t) = 
b a (t)Ji°\ First, solve for dynamics of J[°\ as follows 



tfi(f) 



exp 



ij[ 0) J dt'b l (t') 



ing fields" for 1 < a < d , Ea>i 6 °(^i(*) J « t/i(t) 



Determine new "driv- 
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J2 a >i b a (t)Ja°K Note that no J\ u> -dependent terms are form 
produced in the solvable case. Repeat the procedure by 
solving for dynamics of in the new field, C^(t) = 

, etc., until the last generator is 



?(o) 



exp 



-i4 0) j dt'b 2 (t') 



eliminated and the problem is solved in the following fac- 
torizcd form 



U Ao (t) = f[U a (t). 



(A2) 



0=1 



The difference between this factorized solution |B01 IBT] 
obtained from the standard Schrodinger equation ( 10 ) 



and a solution to the DSBE ( 15 1 is that the latter would 



determine the evolution operator as a single exponen- 
tial directly, effe ctiv ely implementing the BCH identity 



in the product ( A2 1 , U^ (t) = exp 



However, the solvable structure remains in the DSBE as 
well, which should factorize into a chain of ordinary dif- 
ferential equations. An example of such a solution for 
the nilpotent algebra, ()3, is presented in Sec. IV A 



3. Simple-algebra components 

The standard Lie-algebraic terminology becomes a bit 
misleading, when we consider the simple algebra compo- 
nents, which in fact are the most complicated, because a 
simple algebra is an algebra whose generators in a sense 
are all "tied together" via commutators (any element in 
a simple algebra is expressible as a commutator of two 
other elements). In this case, the factorization trick is 
not useful and the corresponding DSBE has the form 
of coupled matrix differential equations, as opposed to 
a chain of ordinary differential equations for solvable or 
nilpotent algebras. On the other hand, the simple al- 
gebras have a very rigid structure. A rank-r simple al- 
gebra can always be represented as an Abelian subalge- 
bra of r-mutually commuting generators, hk [analogous 
to J z in su(2)] and "raising/lowering" operators E a and 
E-oc [analogous to J in su(2)] labelled by allowed root- 
vectors in the r-dimensional Euclidean root space and 
such that [hi, E± a ] — ±ctiE a . The possible root vectors 
in the simple algebra are strictly constrained and so are 
the commutation relations, [E a , £73], which gives rise 
to a complete classification of the simple algebras based 
on the analysis of root spaces and Dynkin diagrams, as 
discussed in the vast literature on the subject (see, e.g., 
Ref. |BS]). 

If the dynamics of the simple component (s) of A have 
been determined, to treat the remaining non- universal 
solvable component, Aq, if any, becomes a straightfor- 
ward exercise that can be accomplished via the same 
factorization trick. Indeed, assume that the mutually- 
commuting evolution operators are known explicitly, 
U Av (t) = exp • JO)]. Then, the total evo- 

lution operator can be written in the following factorized 



U A (t) = 



u Ao (t), 



(A3) 



where the last factor, U Ao (t) = exp [— i<&o(t) ■ J^] , sat- 
isfies the Schrodinger equation in "interaction represen- 
tation," id t U Ao (t) = b (t) • J (0) UAo(t), witn the rotated 
dynamic fields, determined via 



b a (t) 



it 



Ji 0) 



Note that since [A„, Aq] C Aq, the "interaction repre- 
sentation form" of the Schrodinger equation closes within 
the solvable component and the remaining exact solution 
in Aq can be obtained from a chain of ordinary first-order 
differential equations by a variety of standard methods. 

Therefore, classifying possible quantum dynamical sys- 
tems and the corresponding DSBE reduces to an anal- 
ysis of the corresponding equations within each of the 
classical algebras [i.e., A t ~ su(l + 1), Bi ~ so(2l + 1), 
Di ~ so(2l), and C/ ~ sp(l) ] and the five exceptional 
algebras, G2, F4, E 6 , E?, and E$. These possibili- 
ties exhaust completely all universal quantum dynamics 
that can possibly arise from an arbitrary non-equilibrium 
Hamiltonian matrix (modulo a non-universal solvable or 
nilpotent part, whose solution is comparatively trivial as 
discussed above). Apart from fundamental mathematical 
restrictions on the structure of the simple components, 
physics provides further constraints: if a problem can 
be formulated in a physical L-dimensional Hilbert space, 
then we must demand that the group, G A , be a subgroup 
of SU(L). Eventhough, as we discussed above, quantum 
dynamics is completely decoupled from the Hilbert space, 
the very existence of it and the unitarity requirement pro- 
vide certain constraints on types of allowed evolutions, 
U(t) C G A C SU(L m i n ) (here, L m i n can be taken to be 
the minimal dimension of a faithful unitary representa- 
tion). 



4. Generalizing theory to certain 
infinite-dimensional representations 

Our original discussion has descended from the text- 
book Schrodinger equation ([I]) with the assumption that 
the Hamiltonian be a finite-dimensional Hermitian ma- 
trix. This led to a finite-dimensional Lie algebra, which 
by construction had a finite-dimensional faithful repre- 
sentation. This in turn led to the conjecture about the 
global exponentiation of the algebra onto a compact Lie 
group. For our purposes motivated by applications to 
physics, it is the existence of a global exponential map 
that is crucial, while the compactness of the group is 
not germane and in fact there are many theories involv- 
ing infinite-dimensional representations where it is man- 
ifestly not so (e.g., the harmonic oscillator). 
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The global covering problem in its full generality is 
an old and very complicated mathematical question |92j 
and we are in no position to comment on it here. How- 
ever, we outline a class of models, where global expo- 
nentiation should hold irrespectively of the compact- 
ness of the dynamic Lie group, Ga- Consider, a Lie 
algebra, A = span { J 1; . . . , J^} with structure con- 
stants, f ab c , which admits global exponentiation. Form 
parameter-dependent linear combinations of generators, 

d 

J a = C a h {e)Jb- This leads to parameter-dependent 

6=0 

structure constants, / a ;, c (e). Taking the singular limit, 
e — > 0, yields in certain cases a new algebra A = lim A(e) 

£— >0 

distinct from the original algebra. This procedure is 
called a Lie-algebraic contraction (a famous example of 
it is provided by the Galilean group viewed as a non- 
relativistic limit/contraction of the Lorentz group as the 
speed of light is taken to infinity, e = 1/c — > 0). An- 
other example more relevant to the ongoing discussion 
is provided in Sec. |IV C[ which briefly reviews how the 
harmonic oscillator algebra, f)4, can be obtained via con- 
traction from the su(2) spin algebra plus an Abelian 
u(l)-component. Both u(2) and t)4 allow global expo- 
nentiation, eventhough the latter exponentiates in a non- 
compact group with no finite-dimensional unitary repre- 
sentation. We want to view this fact as a manifestation 
of the general conjecture that if an algebra, A can be 
contracted from a Lie algebra, A = lim A(e) that allows 

e->-0 

global exponentiation onto a dynamic Lie group, then the 
former also does. In the context of models of relevance to 
condensed matter, this conjecture implies, in particular, 
that a variety of lattice boson models would allow global 
exponentiation even in the absence of a finite-dimensional 
matrix unitary representation. 



sentation is to express the evolution operator as follows 



Appendix B: Classical Manifold and Path Integral 
for Quantum Dynamics 



The interesting question reviewed in this appendix is 
how to generally obtain the classical manifold, M, where 
classical motion is taking place. This is best elucidated 
in the context of the path integral formalism that can 
be generalized to describe dynamics descending from a 
generic Lie algebra and making use of generalized co- 
herent states [HI |93] introduced by Perelomov [4] and 
Gilmore [5] for an arbitrary Lie group. 

The usual path-integral construction [94] calculates av- 
erages or transition amplitudes directly and therefore 
must include information about the Hilbert space. This 
is despite the fact that the "bulk" of the path integral cal- 
culates the unitary evolution, which, as we just saw, is de- 
coupled from both the initial conditions and the Hilbert 
space. The first step in building a path-integral repre- 



U(t) 



— i f dsb a (s)J a 

Te ° 



lim 



N- 



N 

n< 

n=0 



( B1 ) 

where t n = tn/N represent time-slices. Then a particu- 
lar representation is chosen, Ti,[y4], and a representation 
of unity in a proper basis, e.g., the overcomplete basis 
of coherent states, is inserted between each pair of expo- 
nentials in the product (Bl ) for the given representation. 



This procedure is based on the existence of a set of states, 
\g) G Hil(L), parametrized by the group, g G Ga, such 
that 



h= / dn(g)\g)(g\, 



(B2) 



where dfi(g) is a de Haar measure in the group, Ga- Be- 
ing parametrized by the continuous Lie group, the set of 
states, \g), is way overcomplete. However, as emphasized 
by Perelemov [4] , a part of the group integration is always 
redundant and in fact a more convenient parametrization 
exists, still overcomplete, but where the basis runs over a 
smaller manifold, Ai, to be identified with the manifold, 
where classical dynamics is taking place for the general- 
ized Bloch equations (53). 

The Perelomov construction is as follows: take an ar- 
bitrary normalized state, IV'o); m the Hibert space. In 
the context of the dynamical quantum system at hand 
0, the initial condition provides the most natural choice 
of such a state, i.e., |V>o) = |^(0)). Define a subgroup 
H ( | -00 ) ) C Ga as a set of group elements, h G Ga, such 
that 



&|Vo) 



i4>(h) 



l^0> 



(B3) 



i.e., whose action on the initial condition reduces to a 



multiplication by a complex phase factor [in Eq. ( B3 1 , 
h = TlQi)]. This maximum stability subgroup, H (|V>o))i 
defines a coset M. = Ga/H, which does not necessarily 
have a group structure, but retains a topological struc- 
ture. The coherent states of Perelomov, |M), are an over- 
complete basis in the Hilbert space parametrized over the 
manifold, MeM = G A /H. 

As to the path integral itself, its precise form depends 
on what we wish to calculate using it. A rather generic 
structure is the transition amplitude between two coher- 
ent states, which we can call a Green function, 



M(t)=M 2 



G(M 1 ,M 2 ;t) 



[VM(t)]exp{-iS[M(t)]}. 



M(0)=Mi 

(B4) 

Here M(i) C M C R d ~ A*, i.e., the generalized co- 
herent states lie in the manifold M. embedded without 
crossings into the d-dimensional Euclidean space, K d , as- 
sociated with the underlying algebra. Since dim (Ga) — 
dim (A) = d and A4 = Ga/H, the dimensionality of 
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M. is always smaller than d. Quite generally, we expect 
dim (M) < dim (.4) - IHcmk (A). E.g., the Bloch sphere, 
M = S 2 = SU{2)/U(1) ~ S 3 /S 1 is two dimensional and 
is embedded in the three-dimensional space associated 
with the su(2) algebra. 

If we follow the usual logic in deriving the path 
integral and postulate on physical grounds the con- 
tinuity of the corresponding trajectories on A4, 
the explicit form of the action will take the form, 

S = S H [M(t)] + Skerry [M(t)], where the first Hamilto- 

t 

nian term i sS% [M(t)] oc J dsh(s) ■ M(s), and the second 

o 

t 

Berry phase term, S Be rry[M(i)] oc f ds{M(s)\d a \M(s)), 



incorporates commutation relations and/or topolog- 
ical effects such as the Wess-Zumino term in the 
spin path-integral. The exact form of the Berry 
phase action is determined by the structure con- 
stants of the algebra. E.g., for spin algebra with 



pbc = £ abc^ it ig the Wess-Zumino term, S B crry[M(t)] = 

iS J Q du J dsf abc d u M a (u,s)d s Mi,(u,s)M c (u,s), where 
o 

the dynamic magnetization, M(i), has been assigned an 
auxiliary argument u, in a non-unique way, but such that 
M(0,t) = M and M(l, s) = M(s), therefore describing 
a string sweeping the sphere over time jTQl [SS]. An 
interesting question is to see if the topological terms 
arising from dynamics in other compact groups may 
be cast in such a form, being expressed through the 
structure constants of the algebra and measuring a 
volume covered on the manifold in the course of a 
particular realization of classical evolution. 

Path-integral lore due to Feynman suggests that if we 
now introduce a set of local coordinates, (9\, . . . , 0d), in 
M. [where D = dim M. is the dimensionality of the classi- 
cal manifold] , which is always possible since it is a topo- 
logical space, then the minimization of the action with 
respect to Qi will reproduce the classical equations of mo- 
tion from Sec. [vl i.e., generalized Bloch equations (53 1, 



but now manifestly defined on the classical manifold, M. . 

Finally, let us make the following speculative remarks: 
If we consider all possible Hilbert spaces associated with 
physical representations of a given algebra, A, we shall 
obtain a set of classical manifolds. In complicated higher- 
rank algebras, such as SU(3), they may in principle differ 
in topology and even dimensionality |18j . However, there 
exists perhaps a reasonable "estimate" on their possible 
dimensionality as follows, dim [Gj\/U(r)} < dim(M.) < 
dim [G A /U r {l)\ =d-r, where r = SKank (A). All pos- 
sible classical manifolds, A4, can be embedded into the 
Euclidean space, M. d ~ A*, which is the complete phase 
space of the classical problem ( 53 ) for the averages. If 



we parametrize all manifolds in some coordinate system 
in R d and use certain units along the axes, which nat- 
urally should be related to the Planck constant, ft, we 
shall obtain a variety of possible classical manifolds that 
do not fill the entire phase space K d , but select the re- 
gions of dynamics for the averages, which are allowed 
by quantum mechanics. These manifolds also should re- 
spect and actually follow from the constraints imposed 
by r independent Casimir invariants defined by the al- 
gebra. E.g., for su(2), the set of allowed dynamic man- 
ifolds are the concentric Bloch spheres in R 3 with the 

radii, |M| = |, ft, This is a way to visualize the 

quantization constraints imposed by quantum mechanics. 
On the other hand the Bloch equations (531 do not in- 



clude the Planck constant and are formally not restricted 
by any quantization constraints. We can take arbitrary 
initial conditions for M(0) in M. d , and even if they con- 
tradict the quantization constraints, the classical Bloch 
equations would remain solvable. The "paradox" may be 
resolved by taking the quasiclassical limit in the quantum 
problem, ft — > 0, which in the cartoon geometric picture 
proposed above would imply that the set of allowed man- 
ifolds fills up the entire phase space in this limit. E.g., 
the spacing between the Bloch spheres corresponding to 
different representations of the spin would disappear in 
the quasiclassical limit. 
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